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ABSTRACT 


An automated cloud analysis program was developed and 
established on the SPADS computer system at the Naval Envi- 
ronmental Prediction Research Facility (NEPRF). The pro- 
gram evaluates GOES visual and infrared satellite imagery 
simultaneously. The analysis method produces information on 
cloud types, cloud amount, precipitation intensity, and 
cloud *op height and temperature through use of threshold 
tests of radiance, texture, and temperature. A review of 
current work on the evaluation of satellite information by 
computer and by manual analysis is included. 

A maritime region 460 X 462 nautical miles in size was 
selected for test analysis. Th2 satellite imagery was manu- 
ally evaluated and ccmpared to the computer generated out- 
Pers. Reasonably good patterns 2f cloud types, precipitation 
ama. cloud amount were produced by the computer, although 
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I. INTRODUCTION 


The personnel operating Naval ships, aircraft and weapon 
systems need near real-time subsynoptic and mesoscale mari- 
tine weather information to operate effectively and safely. 
Current conventional cbservational data (from surface , raw- 
insonde, and aircraft observations) do not adequately cover 
the maritime region cr provide for mesoscal2 analysis of 
weather phenomena. Specifically reeded are detection and 
monitoring of any adverse weather conditions. This need 
pertains to phenomena ranging from mesoscale (10-100 km) to 
synoptic scale (100-1000 km) which are poorly resolved by 
conventional observations. Satellite observations can solve 
this problem through the acquisition of global high density 
data which are aveżlabie regardiess of local influences of 
geography, surface conditions, or local politics. 

Another problem is the need to provide an analyzed prod- 
Met tO the user in a short (1-3 hour) time frame. Manual 
analysis of satellite data is time consuming and tedious, 
and it does not provide the information in a real-time 
fashion. An interactive comput2r system, on the other hand, 
ean collect, analyze, and produce the required products 


N 


thin the required time limit. 


| ۸۰ 


N 
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This research effort concentrated on T development of 
techniques to specifycloud tyoss, cioud amount, cloud 
height, and precipitation intensity with an emphasis on the 
maritime regions from satellite images. Information on 
critical weatner parameters such as low visibilities, ceil- 
ings, precipitation presence, and intensity can be derived 
from these cioud characteristics. This work uses Geosta- 
tionary Operationai Environmental Satellite (GOES; Visual- 
Infrared Spin Scanned Radiometer  (VISSR) digital satellite 
data from visual and infrared channels which have a one-half 
hour temporal and 0.5 to U nautical miles (n ui.) spatial 
resolution. The data were received ani processed at the 
Naval Environmental Prediction Research racility (NEPRF) in 
Monterey, California using the Satellite Data Processing and 
Display System  (SPADS). This interactiva computer system 
was developed at NEPEF for uss by the operaticnal Navy to 
give real-time information about weather ohenomcna. One sys- 
tem is now operational at che Naval Zastern Oceanography 
Center (NEOC) at Noriclk, Virginia. 

The objective of this thesis is to astabiish and evalu- 
ate a computer program on the SPADS to analyze GOES iméges 
to derive cloud types, precipitation intensity, elend 


height, and cloud amounts. 125 22217171517 sons Soo an zn-depen 
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study of previous research and analysis techniques was 
mnitiated, and, pesg this study, useful algorithms and 
techniques were gleaned. These algorithms and techniques 
were combined to produce the ter computer program.  Pre- 
liminary tests of the computer progran were completeä on tne 
SPADS system at NEPR?. 


Chapter II begins with ar 


(D 


view of previous studies. 
This includes áiscussions of four techniques of cloud type 
classification; (1) cloud index (Harris and Barrett,1975, 
1978), (2) spectral analysis (Liljas,1981a,  1981b), (3) 
3-dimensional nephanalysis (Fy2, 1978), and (4) two-dinen- 
Sional histcaram (Plat*,1981). A discussion of the techni- 
ques for determining cloud height and cloud amounts used by 
Reynolds and Vonder Haar (1977) and othars, and of proce- 
dures far identifying precipitation areas and their associ- 
ated intensities is also covered. Chapter III describes the 
set of algorithms selected for implementation op SPADS. 
Chapter IV presents an outline of the test procedures used 
and reports the results of these tests (including detected 
problem areas). Charter ۷ summarizes the research efforts 
and problems and provides suggestions for future changes and 


additional research. 
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A. INTRODUCTION 

Satellite visual (VIS) and infrared (IR) inagery can 
provide a wide variety of weather information, including 
estimates of cloud type, cloud height, cloud amount, and 
precipitation intensity. In fact, a satellite nephanalysis 
can produce a spatial summary of these important cloud fea- 
tures. An analyst examines the visual brightness to derive 
estimates of thickness of the cloud and the cloud texture to 
differentiate between cumulus and stratus type clouds. By 
using the IR imagery, cloud top heights can be estimated; 
the brighter the cloud image, the colder; therefore, the 


higher it is. With combined visual and infrared, one has a 


greater ability tO Prek 011 hs -Cloud types and, 


furthermore, to estimate the probapility of precipitation. 
One of the main problems 2ncountered in the 1960's and 
early 1970's was the lack of standardization of rules gov- 
erning the nepnanalysis construction. Bach country, <1] 
sometimes each analyst, has used subiective methods for pro- 
ducing cloud analyses. Some of the specific prchiems aris- 


ing from lack of standardization, as discussed in Harris and 
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Barrett (1975), have been: no minimum size is established as 


reference for inclusicn or exclusion from analysis; cirrus 
clouds are rarely identified in niddla latitudes; 
significant cloud areas are subjectively delineated at tne 
discretion of the analyst; only four cloud categories are 


assigned (which consisted of unequal percents of areas); anċ 


no standard nephanalysis construction rules have been 
published. 
With these problems in mind, Harris and Barrett (1975, 


1978), Barrett and Martin (1981), Platt (1981), Fye and 
Logan (1977) and Lil jas (1981a) „, to name a few, have been 
working toward the goal of producing a recommended standard 
for cbjective nephanalysis of sateliite VIS and IR imagery 


through manual and computer analysis. The following sections 


review these representative research efforts in tiis 
endeavor. The discussions are divided into three areas of 
Eoncer» to this thes is; cioud typing, cloud neight and 


amounts, and precipitation. 


B. CLOUD TYPING 


Basically four distinct approaches tc effective cloud 


t1 
(D 


(1) 


typing have been proposed in the literatura. They a 


Mound index (Harris and Barrett, T975; 19/8: Barrett and 
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Harris,1977), (2) the U.S. Air Force 3DNEPH program (Fye and 
Logan,1977; Fye, 1978), (3) Spectral  cross-correlation 
(Liljas, 198la, 1981b), and (U) two-dimensional histogram 
(Platt,1931). These methods are reviewed below. 
1. Cloud Index 

Harris and Barrett (1975) recognized the need for 
Standardization and fcr increased information content. Asa 
result, they developed a cloud indexing system which pro- 
vided specific guidarce and rules for analyzing satellite 
visual imagery manually. Por clarity and understanding, 
they proposed tne analysis be composed of three separate 
layers of analyses; two descriptive and one interpretive. 
Fig. 1 is the flow chart of the stages in the construction 
of this nephanalysis. The authors developed a key (Fig. 2) 
which includes details on the following cloud features; (1) 
percent of cloud, (2) cloud type ( e.g. cumulus), (3) cloud 
structure (including size, shape, and pattern), and (4) 
Meaqnt of cloud tops. The key also gives interpretive 
weather phenomena. 

The first analysis layer contains information on 
geography and cloud cover. The cloud cover is divided into 
five equal percent categories with a minimum size reported 


of a 2 1/2 degree square. 
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The second descriptive layer includes cloud type 
and cloud structure. They are classified in the following 
order: major cloud systems, definite boundaries and finally 
indefinite boundaries. The analyst usss cloud brightness 
and texture to discern six types of ciouds (see Fig. 2). 

The final analysis layer produced is the interpreta- 
meer OL tne cloud features using the information from the 
first two layers. The types of cloud features considered 
important are dispiayed in Fig. 2. All three analysis layers 
combined were designed to be completed within one hour by a 
skilled analyst. 

Barrett and Harris (1977) extended this procedure to 
include infrared imagery for the purpose of supplementing 
the visual analysis. This added capabilties for nighttime 
and hich latitude analyses. They mention three possibilities 
for analyzing IR imagery; (1) use longstanding internation- 
ally accepted codes and symbols, (2) follow the visual pro- 
cedures previously ovtlined, or (3) use a new procedure 
Which recognizes and represents the special properties of 
the IR images. The first two methods would be directiy con- 
parable to the visual analysis. The third approach would 


meer 3135353160531 information to allow identification of 
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Satelite image 


Outline and denote areas of 
unform % cloud cover in the order :- 
a) Major cloud systems 
b) Cloud areas with definite boundanes 
c) Cloud areas with indefinite boundaries 















Subdivde cloud cover areas on the 
basis of specified cloud tvpes.: 
following the same order of analysis 












Identify and portray clearly— 
defined pattern and structural 
features of the cloud fields 


Review, and revise 
where necessary 









Depict recognizable synoptic 
weather systems and subsystems 
in the lower troposphere 


Depict recognizable synoptic 


weather systems and subsystems 
in the upper troposphere 


Complete the analysis with additional 
descriptive terms where appropnate 
Satellite nephanalysıs 





Figure 1. Manual Nephapalysis Flow Chart (from Harris < 
Barrett, 930r -1 s 
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Figure 2. VIS Key to be Used in Manual Nsphanalysis (£rom 
Harris and Barrett, 1975) 
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cloud height through cloud brightness svaluation and a 
Bpecitic idenċificaticn of cloud textura. 

The IR analysis technique identifies the following 
cioud features (see Fig. 3): (U perzent of cloud cover 
(same as VIS technique), (2) brightness of cloud, (3) cloud 
mexture, and (4) size of clouds (same as VIS), plus inte-- 
pretive categories with eightrsen of then the same as VIS 
(heights and three cloud descriptions are not included). 
eer ett and Harris (1977) felt nme s of the three methods 
could produce an anelvsis in 2 single hour if done by a 


trained analyst. 
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Figure 3. 


AT 


e Used in Hanual Nephanalysis (from 
Harris, 1977). 


After establishing a manual mechod for both VIS and 
IR satellite nephanalysis, Harris and Barrett (1978) devel- 
oped an automated objective nephanalysis through a cloud 
recognition scheme. They used an approach which examined 
quantitative cloud brightness and texture within a small 
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a of the whole picture. ihe. 20559 Char: in Fig. 4 out- 
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es their approach. Cloud brightness was assessed in each 
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subarray by summing the brightness counts of only cloud pix- 
els and then dividing by the number of cloud pixels (Ea. 1 
and 2, Table. I). Cloud amount was calculated in each subar- 
ray by summing the number of cloud pixels and then dividing 
by the total number of pixels per subarray. This number was 
then multiplied by 100 to obtain a percent of cloud cover 
qe 3 ire Table 1). 

Their assessments of? cloud texture used statistical, 
instead of structural, measures of local variation in image 
density. One measure they employed was the standard devia- 
tion of each subarray (5x5 pixel array) of density vaiues, 
as calculated with Eq. 4 (Table I). 

Their other criterion was wector dispersion, which 
required that he density values be ‘treated as a set of 


adjecent triangular planes. The dispersion in three dinen- 


Sions of the normal tc these planes gave another measure of 
texture. Pag. 5 sh¥ews a graphical representation of this 
Method. Eq. 5, 6, and 7 in Table I wers used to calculate 


these estimates of texture. 
Once all tnese vaiues were caiculated, boundaries 
between cloud types had to be established. Three categories 


^ = 


eet or layered clouds  (stracirorm), 


Y 


were selected; (1) 





(2) cellular/tower clouds (cumulus) and (3) broken/nottled 
cloud (stratocumulus or mixed) (Harris and Barrett,1978). A 
three-dimensional decision spaca was sət up using two tex- 
ture measurements and one brightness value. Those areas that 
had a brightness above the threshold were included in fur- 
ther analyses; those below wer2 not. Fig. 6 shows the deci- 
sion space based on the two texture msasurements. The 
discriminant lines were developed through analysis or a 
training set of Defense Military Satellites Program (DMSP) 
sample brightness and texture measurements of each cloud 
category. 

In assessing the accuracy cf this method, Harris and 
Barret+ (1978) found several problem ar2as. For sxampie, 
towering cumulus in post frontal areas (esvecially) were not 
identified because of the field of view (FOV) size and snow 
areas were included as clouds. But the overall accuracy was 
meet good wath a greater than 72% correct classification. 
Harris and Barrett (1978) feit that several procedures would 
improve the output of the model. Among those were increasing 
e resolution of the data by decreasing the size of scan- 
feng spot, decreasing the subarray size, normalizing the 


image brightness for changes in sun viewing angle, and 
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increasing the texture parameters to permit identification 


of more cloud classes. 
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Figure 4. Flow Chart for Automated Objective Nephanalysis 
E (from Barrett and Harris, 1978) p 5 


Dinensicnal Nephanalysis (3DNEDA) 


The United States (U.S.) Air Force developed an 
automated three-dimensional cloud analysis called the 


EEDNEPH model", It uses satellite imagery from the DMSP and 


TIROS-N satellite series together with conventional in situ 
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a) 


b) 





Figure 5. Graph of Veczor Dispersior Techni que. (a) 
S resents a smooth surface (Stratiform clouds) 
l- represents a rouqh surface (cumulus 
clou P Normal wvectasrs ror each case are 
included to the E of the diagzams (from Harris 
and Berrett, 1978). 
data from aircraft, surface, and upper air cbservations. A 
cloud analysis is procuced in a real-time fashion regardless 
of ice, snow, desert, or other geographical anomalies (Fye 
Eud Logan, 1977). The satelliras provids high density data 
over the entire world, thus permitting a global analysis of 
Metas. The horizontal resolution of the output is 25 n mi. 
The data array has a pixel resolution of 3 n mi. 
The satellite data are first formatted, rectified, mapped, 


and stored on a computer storage device before running the 


3DNEPH progran. ۸1 :غه د‎ 1١6و‎ ٢1 98 COrrects ror the effects 
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pha ae (from Harris and 
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of limb darkening before storage. Linb darkening is the 
effect of a decrease in brightness as the edges or "limbs" 
of the earth scan are approached. The data are not normal- 
ized with respect to light (VIS) or longwave radiation (IR) 
(Fye and Logan, 1977). Each visual and infrared pixel is 
represented by a single grayshade value in the range fron 1 
8:963. in the visual, 1 is dark and 63 is white. Infrared 
pixels represent temperature in degrees Kelvin (K) from 210 
K for a value of 1 to 310 K for 63; sach grayshade repre- 
Sents a 1.6 deqree change (Fye,1978). Tke preliminary pro- 
cesscr calculates the average grayshade value per 25 n mi. 
square using Eq. 1 (Table II) and then uses this to get the 
Eoasbil2ity within the square using E3. 2 {Table II) 
(Fye,1978). Both these calculations are performed sepa- 
rately for IR and VIS images. The VIS and IR grayshade val- 
ues and variabilities are used to identify the various cloud 
types. Fig. 7 and 8 represent the cloud typing algorithms 
used. They were constructed empirically by comparing nany 
grayshades and variabilities to corresponding satellite 
imagery and surface data (Fy=2,1978). The operational algor- 
ithns were derived from the figure and then modified to 


Bbr32n optimum results. One obvious problem with electing 
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to use the model {given the available documentation) is the 
lack of graduated values on the axes of these two figures. 
Only relative values of warm/cold or light/dark are included 
in the report. Also lacking in the Fye(1978) report are any 
details of the current algorithms used for cloud typing. 


TABLE II 


3DNEPH Statistical Equations for the Average Grayshade and 
Variance (from Pye, 1978) 


N G = average gray shade 
(1) G = 1/N x G. , G>0, N>0 - - individual 3 nmi area 
T : grid gray shade 
N = total number cf G, per 
25 nmi grid space 
N 
(2) V = 1/8 Y le, -&| ,c>0, N>0 V = variability within the 
ml 25 nmi grid space 





Mie third Classification Wethod is the spectral 
cross-correlation method developed by Erik Liljas (1981a) 
for Use On an interactive computer. Liljas used TIROS-N and 
National Oceanic and Atmospheric Administraion (NOAA) 6 
satellite data from three spectral channeis; channel 1 
(VIS), channel 2 (near-IR) and channel 4% (IR). During the 


summer time, Gifierent cloud types, land and water surfaces 
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Figure 7. The Decisicn Space of the 3DNEPH Program. IR and 
and VIS grayshades are used to determine the 
cloud type (from Fys, 1978). 
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were analyzed to detect their unique spectral signatures. 
The results were used to build a three dimensional 
paralielepipedon classification, Fig. 9. 

Two satellite types were used for the test of the 
model; polar orbiting (TIROS-N and NOAA-6) and geostationary 
(Meteosat). The satellite data used for the test of the 
model were geometrically corrected, but they did not need to 
be ncrmalized because the mid-day pass of “he satellite was 
used (approximate sun elevation SI 45 degrees). 
Normal!zation corrects for the radiance variation caused bv 
the sun location (elevation). Also the amount of data was 
reduced from 1,0240 digital levels to 256 levels, which was 
considered adequate resolution for the study. Data from 
channel 2 were primarily used to separate land and water 


areas. Channels 1 and 4 were the principal cloud =ype 


Comparison of the spectral cross-correlation method 
With synoptic observations showed that the computer method 
gave correct indication of cloud *ypes, and also an improved 
cloua division. As shown in Fig. 10 (see Table III for sym- 
bol definitions of cloud types), several cloud tyves are in 


the same boxes; Liljas recommends the use of texture 
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observation to discern between these cloud types. This was 
done manually and entered interactively. Liljas (1981b) also 
reports that good infcrmation on cloud structure of weather 
systems can be retrieved. He recommends the use of nultiple 
parallelepipedon classifications, each classification asso- 
ciated with a different sun elevation instead ot 


normalization of the reflected imagery. 
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Main and Cloud Types: 


1. 


Cumulonimbus 


Nimbostratus 


Cirrostratus 


Cumulus congestus 


Stratocumulus 


Haze/Stratus 


Land 


Water 


TABLE III 
Cloud Classification to SE AE Pig- 10 (from Liljas, 
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Storm cloud with high top 
Squall cloud with scattered showers 


Large vertical thickness 
Rather low topside 


Dense cirrostratus 

Cirrus 

Thin cirrus over water 
Dense altostratus 

Thin altostratus over water 
Thin altostratus 


Dense altocumulus 
Large piled up cumulus 
Rather small and flat cumulus 


Dense stratocumulus 
Ordinary 

Slightly piled up cumulus 
with clear areas in between 
Very dense haze/stratus 
Dense haze/stratus 


Ordinary haze/stratus 
Cumulis humilis 


Haze over water 
Planting season spring or autumn 
Warm green season 


Cold 
Warm 
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Another spectral method was developed by Piatt 
(1981) using a two-dimensional histoqram representation of 
cloud types. He selected satellite data with typical cloud 
fields to develop a histogram decision space representation 
of cloud types. Platt concentrated on evaluating cloud sys- 
tems instead of isolated clouds and therefore his results 
are biased toward syncptic features. 

The measured visual radiances were compared with the 
direct radiance from the solar disk on the satellite to pro- 
duce a caiculated isotropic albedo per 2.5 km. These aibe- 
does are expressed as a fraction from 0 to 1. The infrared 


Bad-rances are corverted to blackbody brightness t 
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MEA» from 193.3 to 301.5 4. A histogram was then 


Ui 
Ó 
rh 


constructed by sorting these data into 983 different bin 
temperature and albedo. Each area evaluated was 2pproxi- 
mately 500 X 500 km, with roughly 40,000 values per histo- 
gran. The total Den counts wera plotted on a 
two-dimensional grid and then contoured by same bin count 
(see Fig. 11). This cloud histogram representation could 
then be used to classify cloud types by matching plots of 


the satellite data. After classifying various types of 
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clouds, Fig. 12 was compiled to give an idealized depiction 
of the location of clcud types by albedo and temperature. 
Platt (1981) discussed several factors that tend to 
reduce the usefulness of Fig. 12. One cloud problem arises 
when there are either breaks in a cloud deck, or partially 


filled field of view (FOV) in the satellite view area. 


(+ 


tionaily, clouds or the sama Optical depth, but exhibit-‏ د 


fu 


ing breaks will have a ifferent aibedo measurement thar an 
unbroken cloud of the same optical depth (stratocumulus in 
Fig. 14 demonstrates this). Another factor that changes the 
albedo and temperature values is variation in cloud top 
height. Ali cf these factors create ambiguity in determining 
the cloud type, amcunt, optical depth, and cloud top 
temperature. As a result, the author recommends +har addi- 


tional work should be done with these results before using 


the two-dimensio ral histcgram approach to identity clouds. 
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C. CLOUD HEIGHT AND AMOUNT 
Three methods of deriving clcud height and amount have 
been proposed in current research. Two methods are an exten- 
sion of those reviewed in the cloud typing Section; cloud 
index {Harris and Barrett, 1975, 1978; Barrett and Harris, 
1977) and the U.S. Air Force 3DNEPH (Fye and Logan, 1977; 
Fye, 1978). In addition, a new bispectral technique 
(Reynolds and Vonder Haar, 1977) is discussed. 
1. Cloud Index 
In two separate studies, Barrot- and Haris (1977; 
Harris and Barrett, 1975) developed a recommended standard 
cloud neight and amount determination procedure for manual 
analysis of satellite visual and infrared images. The vis- 
ual method combined brightness, cloud type, and synoptic 
situation to determine qualitative estimates of cloud height 
mow, middle, or high level). Later, they developed a sys- 
tem for identifying cloud levels through the dearee of cloud 
brightness in the IR image. IR cloud information specifi- 
Ba ly lends itself +o this Ca k Since the rač arti on 
momperature of the cloud is a function of the cloud-top 


EM tudce (Berret* and Harris, 1977). 
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Some problems arose with cirrus and cloud elements 
smaller than the field of view. In this instance, satellite 
measured temperatures can appear warmer than actual cloud 
temperature, because the satellite is also receiving the 
radiation temperature or the surface or lower clouds. Dus 


prebiem can be minimized by using IR and VIS together to 


Maty cirsms clouds and by using higher resolution data , 
S. Jo, date from the Advanced Very High Resolution kadiometer 


(AVHRR) on NOAA satellites. 

Eg 3 is the key for the IR nephanalvsis and 
inciudes the cloud brightness steps and Fig. 13 is an exam- 
ple of how these are applied. This analysis can be inter- 
preted in terms of quelo 6۹2011 Los (GYa) to high 


[ery 5Drign:) Bazre** and  Herrzis,1977). Cloud amount 


re also delineated in Fig. 2 for 


pi 


determination categories 
MES and Fig. LODE R; It is donre in categories of size of 
cloud =iemenz, روو د2‎ 110227 km. dacris and Barrett's 
(1978) objective automated nethod does not specificaily 


0802255 cicud height, put does caiculat2? an average cloud 


w 


ment Sor each subazray size (5x5 pixels) using Eq. د‎ 


(Table I). 
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Figure 13. 
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2. Three-Dimensicnal Nephanalysi 


(3DNEPH) 

The U.S. Air Force 3DNEPH program's main objective 
is Lo produce information on cloud ut and cloud amount. 
It uses a standard vertical grid of 15 layers, which range 
from the surface to 100 millibars (mb). The IR image is 
used to determine height and element cloud amounts, and the 
visual and infrared are used together to determine total 
cloud amount. TALS Sra ON :s then conbined with sur- 
face, aircraft and rawinsonde data to determine the final 
cloud height and cloud amcunt analyses. 

The preliainary processor, as Bee, in the cloud 


type section, prepares the satellite data for height deter- 


mination. ipad sen, 3-5 90Dts, screens, and combines the 


l-9 
taj 


conventional data received from the Air SEES *“ 1 
weather data base to produce a vertical and horizontal 
3DNEPH gridded layered format. 


Next, the satellite data goes to the visual and 


t4 
(D 
Cu 
j- 
JI 


infrared processors. The VIS processor steps are outli: 
Erg. 114 . Baci ۵ل اوح ړل هده اح دده‎ is direstiy cohpared to a 
grayshade (backgroundysurface brightness) specified for each 
25 n ni. Square. This value represents a clear sky eart! 


radiance for each square. The difference between each دم‎ 
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grayshade and the background brightness is used in Fig. 15 
to decide if there is a cloud or not. Th2 curve in Fig. 15 
represents the brightness threshold for clouds; any vari- 
ability greater than this isa cloud area. This curve 
incorporates the background brightness and the varying 
brightness associated with different surfaces ( e.g., 
desert). The total number of cloud decisions divided by the 
total number of available pixels represents the amount of 
cloud cover in the square. 


The infrared processor incorporates a slightly dif- 


tty 


erent anaiysis method. The majer steps of the processor 


in Fig. 16 asa flow chart. The first step is‏ 226106826 ه 


pu 
1 


to construct a frequency distribution or the 64 raw gray- 
Shades for each 25 n mi. area to detect significant nodes 


(maxima) cf 6 or mere qrayshades or 12 or more adjacent 


if 


des ac 


0 


shades. I2 more then 2 modes are identified, the a 
combined on the basis of relative amplitude (strongest nodes 


are retained) and «he differences between modes. if there 


Lis 


Sno a 
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perer- mods, then the processor 3122381158 to a cloud 


ch 


hreshold temperature determined from appropriate surface 
temperatures. Before any cloud decision is made, 2 correc- 


tion is applied to the IR temperatures to compensate for 
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atmospheric absorption and other effects such as zenith and 
local angle which create biases in the IR data. Fig. 17 is 
an example of a curve used to derive bias correction for a 
particular satellite. The actual bias correction is the dif- 
ference between the best fit curve and the diagonal line of 
perfect correiation. Next a representative temperature is 
selected for each mode (usuaily the coidest in the mode} and 
these data are compared to a gridded, conventionaily derived 
Surface temperature data base to make a cloud/no cloud deci- 
sion. Fig. 18 graphically represents this decision process. 


Ene final cloud decision is made by using Fig. 19. THS 


E 


compensates for rreqularities in the surface temperature 


data base which limit the reliability of the cloud/no cioud 
decision. The cloud amount is determined by -he number of 
pixels in the cloudy mode  (Fye and Logan,1977). Adigscon- 


ally, the IR processcr determines the cloud height by com- 
Baring the coldest temperature in the cloudv mode to a 
real-time upper air data base. 

Othez processors include the conventional data pro- 
Gesscr, the data integration processor, and the manual data 
processor. The conventional data processor takes the sur- 


face, aircraft, and upper air reports and produces a 
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horizontaily gridded, eight mesn data base of cloud informa- 
Exon for 15 vertical layers. These layers of data and the 
output of VIS and IR processors are merged by the data inte- 
01352101 processor. The processor also perfcrus meteorclogi- 
cal consistencv checks. The cloud amounts, heights, and 


types from VIS and IR data ars converted to 


0 


pecastc cloqd 
layer amounts and heights. Conflicts are also resolved bv 
this processor by checking age of data, greatest cloud 


Sewer, layer amount, and cloud top height. Next, the 
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satellite and conventional data are compared and integrated 
with a continuity field to complete the cloud analysis 
(Fye,1978). The continuity field is the previous or latest 
3DNEFH analysis which fills in missing data and parameters. 
The analyst next enters  manuaily encoded input to 


alter the 3DNEPH data base as necessary. The manual data 


U 


processor takes ‘this information and incorporates it into 
the cverall data base. This provides an averue for quality 
control, correction of errors and improvement of the overall 


analysis. 


E 


{ng 
lu 


ispectral 


بابر پوت 


y 


bsspeceral meuhed for @determining cloud height and 


amount was developed by Reynolds and Vorder Haar (1977). 


(D 
Cu 
(D 
م‎ 


Visual and intrar data were each anaivzed quantitatively 


ps 
O 
c 
fs 


over an array of scan spots. This analysis yielded c 
amount and cloud top temperature, from which cloud height 
was derived by comparison to a nearby vertical temperature 
sounding. 

The method was based on manipulation and solution of 
imo equations {Eq. 1 and 2, Table IV). These equations are 
basically ccmpriseä of adjustments for the fraction of 


cloud-covered and clear-sky areas in a single field of view 


and for backbody emissivity of the clouds. One was the 
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satellite visual radiance equation (Equation 1, Table IV) 
and the other was the satellite infrared radiance equation 
(Eq. 2, Table IV). The equations related the amount of cloud 
(Acld in tenths) to either cloud albedo (pcld X Hs) in the 
visual or cloud radiances (Ncld) in the infrared channels 
(note that Acid plus the amount of clear area, Aclr in 
tenths, equaled 1.0). The IR radiance equation also 
accounted for the effect of smissivity (ce) of clouds on 
cloud and surface (Nclr) radiances. In both equations, Aclr 
was replaced by 1.0- Acld. The rest of the terms in the 
equations except two (Acld and Ncld) were available. Average 
values 5f the cloud and surface albedoes were known and more 
exact values can be cbtained in practice "from the actual 
array of satellite measurements" (Reynolds and Vonder Haar, 
1977). Measured surface temperatures were converted to IR 
radiances and used for the surface radiance. The enissivity 
was initially assumed to be 0.9 for all clouds. This left 
only the two unknowns to be deternined. EI. 1 (Table IV) was 
ther solved for Acid and Eq. 2 (Table IV) for Ncld radiance 
which resulted in Eq. 3 and 4 (Table IV). The Acld value was 
Aer i več first from the visual sacellate data. Next, 83 


was used with the simultaneously received IR data to obtain 
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the cloud radiance value. The resulting radiance was then 
convertad to cloud top temperature via Plancks’ equation. 
The final step was to use a local upper air sounding and the 
cloud temperature to obtain the cloud top height. 

Reynolds and Vonder Haar (1977) performed an uncer- 


tainty analysis on this technique to determine on what the 


uncertainty (A) of the cloud amount (Acid) and cloud 
temperature (ANcld) measurements jepenä and how Large an 
uncertainty (A) is allowable. it was shown that A Acid 


depends on the m agnitude and precision of Hs, Ms, pcld, and 
pcir, where ۸٥815 and A pcir will be less than 10% of the 
error.  Ncid cepends on the magnitude and uncertainty of Ml, 
Ed, Nclr and €. 

The two quantities 423-13 and ANcld were evaluated 
using typical mid-latitude values obtained from satellites. 
Table V gives the absolute AAcld error for givenA Hs, Ms, 
and pcld. These results show that the "uncertainties mus“ 
stay below the 10% level i£ reascnable accuracy is to be 
maintained" (Reynolás and Vondar Haar, 1977). This accu- 
racy is possible through use of ground and cloud truth neas- 
urements along with improvements in instrument calibration. 


At present there is approximately a 3 degree Celcius (°C) 
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(1) 
(2) 


(3) 


TABLZ IV 
Algorithms for Cloud Top Temperature and 53313 Amount (after 


Reynolds and Vondar Haar, 
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error in satellite derived temperature which is equivalent 
ZO approximately 5% Ain Nclr. Ncld will have a 16°C er- 
ror i£ the Acid has 24% error (Table V). This error can be 
reduced by better Acld determina tion using higher resolution 


data andi through use cf ground truth measurements. 


TABLE V 
imMEWErrors in Cloud une 55 25 ل‎ 8572 5%) and in Cloud 
Temperature Necla = 2 و‎ Reynolds Ma Yonder Haat, 


55 pclr, IA ocir, 155 pclr 


Uncertainty pcld, Ms pcld Ms ٥614, 5 
AAcld (Acld = 0.26) 0. 09 0.14 0.20 
AAcld (Acld = 0.53) 013 0.21 0.30 
AAcld (Acld = 0.83) 07 029 0.40 

SAME, * 10% ML, 
Hels. € 117 E 


24% Ln 469 40% in Acid 


ANeld 16° 307 
AZ 3 kn 4.2 km 


The bispectral method was tested with NOAA satellite 
coverage in an area that had cloud observations and upper 
air dati available. Three locations were chosen: White Sands 


Missile Range (WSMR), Denver, and Oklanoma City. 1ل‎ SITS 





km area was used because it best corresponded to surface and 
upper air data. Average values of pcld, pclr, and Hs were 
assumed initially and then changed as sun angle and 
satellite viewing angles varied (Reynolds and  Vonder Haar, 
1977). The results showed a rms error of 0.2 in cloud 
amount with a slight underestimation bias. The cloud height 
error with cirrus included was 4.6 km with a -4.2 km bias. 
Without cirrus cloud, the rms arror improved to 0.5 km with 
overestimation of 0.27 km. 

5515 2616م‎ with cirrus is linked to the use of 0.9 
for the emissivity of alli clouds. Because of the ice crys- 
tals in cirrus, its emissivity is less than the 0.9 value. 
Two approaches weze considered by Reynolds and Vonder Haar 
STD o correct this problen. 

One approach was to correlate the visible spectral 
albedo of a clcuá to its emissivity through the blackbody 
temperature. Fig. 20 shows the scatter diagram of the tela- 
tionship between the albedo and the blackbody temperature 
and a second order fit which was used +o determine che cioud 
top temperature for a specific albedo. Ince this cloud «op 
temperature is known the emissivity can be calculated 


SING Eq. Y, Table IV (derived from Eq. 2 of Table IV). 
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Fig. 21 is the resulting curve which relates aibedo to emis- 
sivty. The curve is a second order equation of the best fit 
curve measured in the bispectral method along with a slight 
adjustment to pcir and pcid. Dies دل«‎ the vennor Cr the oid 
method in half (Reynolds and Vonder Haar, 1977). 

The second procedurs for determining cirrus cloud 
Mm'eqht Was 20 compare channeis 4 and 8 of the Vertical 
Temperature Profile Radiometer (VTPR) on the NOAA-U 
satellite. This led tc a reduction by one half in the appar- 
ent error of the oid method. Tabie VI lists the results of 
using a cirrus correction and the VTPR approach compared 


Men non-adjustaent for cirrus. 


DA ESSCIPITATION 

Four approaches for estimating precipitaron from 
Satellites are reviewed in this section. They are cloud 
dex (Barrett and Martin, 1981; Barrett, 1981), iife his- 
Bory (Barrett and Martin, 1981: Griffith and Wocdley, 1981), 
Be-spectral (Liljas, 1981a, 1981b; Austin taza Lovejoy, 
moe, ana cloud modeling (Wyli2, 1979; Simpson and Kigger“, 


MJO9). 
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Figure 20. The Relationship 3etween Albeio and Blackbody 
Radiation. This was derived from NOAA sat -eliite 
measurements of cirrus ciouis over the Gulf of 
california (fron Reynolds and Vonder Haar, 
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Figure 21. The Relationship Between Albedo and IR 
Emissivity from Empirical Studies (from Reynolds 
and Vonder Haar, 1977). 


























TABLE VI : 
0 15215021 cf cirrus correction and VEPR results (from 
Reynolds and Vonier Haar, 1977) 
CIRRUS CORRECTION RESULTS 
Date Uncorrected Cloud Height Bi-Spectral Cloud Cloud Height Observed Cloud Top 
Bi-Spectral Cloud Error (Km) Top Temp. with Error (Km) Temp. ( X) 
Top Temp. ( K) Cirrus Correction (xK) 
8 Jan 75 258 *3.2 248 *1.7 233 
8 Jan 75 261 +3.3 242 +1.1 233 
23 Oct 75 256 +2.2 229 - .5 233 
9 Oct 75 250 *1.6 231 - .4 235 
RMS Error = 2.6 km RMS Error = 1.1 km 
Bias = +2.6 km Bias a +.48 km 
YTPR RESULTS 
Date Spot YTPR Determined SR Determined Observed Cloud Error (Km) 
Cloud Height (Km) Cloud Amount Meight (Km) 
3 Jan 75 33H 114W 8.2 .7 8.7 £5 
8 Jan 75 33.24 114.34 8.6 .7 8.7 * .1 
8 Jan 75 33.5N 116.34 9.7 .6 8.7 *3.0 


RMS Error * 1.75 km 
8ías «+1.2 km 
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1. Cloud Index 

Barrett (1981) developed a cloud index method of 
determining rainfall (called the "Bristol Method’). It is a 
manual method which uses observed rain gage data and the 
cloud index system (see Section B) to predict the amount of 
rainfall accumulated over a spacified time period (usually 
12 hours) in arid squares of a selected size (1/69 to 19 
grid squares). Cloud area and altitude ars also taken into 
account (Barrett and Martin, 1981). Fiz. 22 presents the 
flow diagram of a recent version of the Bristol Method with 
two stages, preparatcry and oparational. Bosch IR and YIS 
are used wherever possible to identify cloud types. Next, 
the cloud types of aach grid square are assigned an index 
(see Table VII). The cloud index is generaliy "evaluated 
differently for ranges of chosen cloud types and for differ- 
ent regions" (Barrett and Martin, 1981). Then mean cloud 
indices are calculated along with mean observed rainfall 
amounts for each cloud field unit. Subsequently, regressien 
techniques are used to relate the observed rainfall to the 
cloud type, and these nelp establish rainfall estimates for 
areas of satellite cells where conventional data are sparse 


(see Fig. 23). Tests show that the accuracy of this method 
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is approximately 75% 


for all cases belonging 


to a specific 


category (Barrett and Martin, 1981). 
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TABLE VII 


Ed ne irfall Probabilities and Intensities‏ إا ياك 
fron Sate lite Pictures ( from Barrett and Martin, 1981)‏ 
3 2 1 
Assigned probabili- Assigned intensities‏ 
States of the sky ties of rainfall of rainfall (rela-‏ 
(nephanalysis cloud (relative scale range tive scale range‏ 
categories) 0 - 1.00( 0 - 1.00(‏ 
Cumulonimbus 0.90 0.80‏ 
Stratiform 0.50 0.50‏ 
Cumuliform 0.10 0.20‏ 
Stratocumuliform 0.10 070)‏ 
Cirriform 0.10 0701‏ 


Clear skies == a 
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The Griffith/Woodley (Griffith and Woodley, 1981; 
Griffith et al., 1978) technique for estimating convective 
precipitation was developed using a cloud life history model 
in which satellite imagery is compared to radar echo life 
histery of the cloud area (Griffith and Woodley, 1981; Bar- 
Bett and Martin, 1981). The basis of the technique is the 
MMWMOENESsS that in the tropics rain occurs with cold 
GE) clouds and the intensity of the rain is a function 
of the stage of development of cumulonimbus clouds. in a 
Slay Of cloucs over florida, Griffith and  Woodley (1981) 


1 تت‎ the  cumulcnimbus cioud aree is related to radar 
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Figure 23. Regression Diagram Used to "Float" Rainfall 
iStinates ger Stell ite Slctur as (ironia ost 
and Martin, 1981). 

echo area (see Fig. 2U). Cloud areas (Ac) are evaluated in 


either VIS or IR satellite imagery and then compared to tne 
radar echo area (Ae). Fig. 24 is a plot of the relationship 


between cloud and echo areas (Griffith 


[o 


t alos; UO) as 
empirically Gerived using measurements and compar- -sons mor 
710223 Clouds Viewed by radar and geosynchronous 
-- The tain cloud threshold (IR) Usario. e 


approach was 253 K. 
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This technique is rully automated (Figure 25) ana 

has been tested in the Florida area and in axtra-tropical 
convective areas (Barrett and Martin, 1981). There are four 


major computer modul=< that together produce a rain map. 


The 





COM. are navigatica, cloud isolation and tracking, cloud 
life history compilation and volumetric calculation of rain, 
and rain mapping (Griffith and Woodley, 1981; Barrett and 


Mertin, 1981). 


(b 


This model has been applied on several time scalés 


and in different regicns as shown in Tabie VIII. The error 


fu 
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fs 


ERES 2stics (bias (3), meot Man Slats srror  (S-ms), 


tt 


linear least squares fit (R, slope, and intercept)) DT 
these applications (Table IX) show that on average the model 
results agree fairly well with the grcund truth data 
(Griffith and Woodley, 1981). The EENS is small except for 
two locations, with an overall average of 0.85 millimeters 
per hour. The bias is close to 1.00 which is the value of 
perfect correspondence. In tha Macas, all correla- 
tions (R) are greater than or «equal to 3.78, where 1.0 is 


perfect correlation. 


3. Bispectral 





Lovejoy and Austin (1979) developed a technique in 
which visual and infrared wavelength patterns are used to 
|1 117۶7 clouds with a high prcbability of rain. These 
clouds have large cptical and physical thickness "as 


measured by «he visible albedo in comparison with their 
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Bae GATE (from Barritt and Martin, 1981). 
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Applications, 


AREA 


Global Atmospheric Research Project's 


-TABLE VIII 


he G 
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an 


Atlantic Tropical Experiment (GATE) 


Florida Area Cumulus Experiment 


(FACE) 


Equatorial Pacific Ocean Climate 


Studies (EPOCS) 


High Plains Experiment (HIPLEX) 


Hurricanes 


Flash floods 


LOCATION _N 
GATE 53 
HIPLEX 15 
(dense) 

FACE 11 
HIPLEX 9 
(sparse) 


HURRICANES 3 


FLASH FLOOD 2 


PERFECT 
CORRESPONDENCE 
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kjoodie‏ م 
woodley, 1‏ 


ss (from 


PERIOD 


June 27 to Sept. 20, 1974 


Selected days June to September 
1972 to 1980 


2-1/4 days November, 1978 
30 days November, 1979 


Selected days, May to July 
1976 to 1978 


Selected storms 1969 to present 


Johnstown, PA July 19-28, 1977 


TABLE IX 
from the Appl 
GE EAC d 

B ERMS R 
9.84 8.79 8.87 
8.77 2.81 8.98 
0.88 8.42 8.95 
8.98 8.55 8.78 
1.08 0.26 1.08 
1.39 0.48 0.99 
1.08 0.00 1.00 


ations Listed in Table 
Gatley, 1981) 

SLOPE INTERCEPT 

9.97 -1.47 

0:62 8.97 

1.18 -9.23 

1.81 4. 27 

0.89 0.18 

0.89 23.93 

1.28 1.00 





height, determined ty the intensity of the IR emmissiont 
(Austin and Lovejoy, 1981). Liljas (1981a, 1981b) followed 
this with a simliar type of bispectral identification of 
precipitation area based on VIS and IR imagery. This sub- 
section will first review Austin and Lovejov's (1981) 
technique and then Liljas* (1981a, 1981b) technique. 

Austin and Lovejoy (1979) first addressed the ques- 
tion, "How well can rain areas be delinsated in coincident 
visible and infared satellite images?" They felt that the 
problem was to find boundaries between rain and no rain. 
Estimation of rain amcunt is assumed to depend on two steps, 
the delineation of the rain areas ard =stimation of rain 
rate once the area is known. The authors used three statis- 


[meer to charecterize errors and fron these were able +0 


t3 


zt WSeepein rein 


00 
rn 
j-e 
ti 
Ui 


estimate rain area and rain rate. h 
estimation was to map a frequency matrix for both rain areas 


and non-rain areas (see Pig. 26 and 27). The ratio of zain 


to non-rain was found and plotted (see Fig. 28) along with a 


line of greater then 50% probability of rain. Premo his, a 
rain map (see Fig. 29) was produced. The next step was ٥م‎ 
Ey oOo quantitatively estimate the rain rate. Austin and 


Lovejoy (19651) tried several tschniques to do this but met 
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with little success. They believe this was due to "...the 
visible and IR wavelengths primarily responding to the rela- 
tive abundance of cloud droplets and not to precipitation 
sized particles" (Austin and Lovejoy, 1981). They concluded 
that a good or even passable guantitativ2 estimate of rain- 
Pall tate is not realistic or practical using visual and 
infrared data. The authors f21t their study supported the 
hypothesis that visual and infrared satellite images are 
good for estimating rain areas but not for estimating 
quantitative rain rates. 

Liljas' .(1981!a) main objective was to appiy an 
apprcach Similiar to that used by Lovejov and Austin (1979) 
for identifying precipitation areas and intensity using 
TIXOS-N satellite imagery, but without using radar to cali- 
Beate it. He decided not to try quantifying intensity, but 
to instead define it qualitatively by categories of light to 
heavy rain. This fcllows thea work of Muench and Keegan 
mae o) on the relation of cloud reflectance and intensity of 
rain. Pursuing their concepts, Liljas (1981a, 1981b) devel- 
oped a relative ده و و‎ picatiom intensity identification 
technique which relies on his cloud typing model (see Sec- 
pon B, Chapter II). CES 52 2221312-84 quickiy to both 


frontal and convective type precipitation and gives the 
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information in a real-time fashion. The first step is to 
identify the cloud types (see Fig. 9). Next, the area on the 
graph representing cunulonimbus and ninbostratus is sec- 
tioned (see Fig. 10) based on the premise that cold, high 
and dense clouds give the highest precipitation rate. An 
additional requirement for the cloud top temperature to be 


ljas, 1981a) and sup- 


colder than -220C was established (L 
ported by investigaticns cf Muench and Keegan (1979). Table 
X gives the thresholds determined by summing the infrared 
and visual digital values that correspond to the precipita- 


tion intensity levels. The method gavə Jood information on 


the distribution of precipitation areas and the relative 
intensity within the areas for both fronta ard scattered 
Showers. In addition, the method detected "intense precipi- 


55531623 in mesoscale cleud clusters between synoptic! 
observations (Liljas, 1981a). 
4. Cloud Model 


کح e yes SS a‏ اص پا 


Wylie (1979) compared <ropical and mid-latitude 


Ft 


cloud area and rain rates to see i the GriffithsWocdley 
Ben Erich et al., 1978) and other techniques could bs 


Bea tO non-tropical cloud features, instead of the trop- 


ical cloud systems for which they were developed. Six days 
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TABLE X 


Threshold Vaiues Describing ee 


ita nsity Levels 
as Applied in Fig. 10 (from L 1a) 


Chl + Ch4 291-310 light rain 
311-330 
331-350 
351-370 
371-390 


2990 very strong rain 


MEMEWO different types cf rain clouds (frontal convective 
End frontal Lco n the Montreal, Quebec, Carada area 
were chosen for comperison to the 1979 GATE data set of 
Out et al. (1979) (see Table TII). 

The rain rates of the two areas are shown ir Fig. 
BO. The GATE rain rates are higher than the rain rates of 
Montreal, even at the colder threshoids. mne Montsesel rain 
rates tend to vary ccnsiderably, especially the convective 
cases (June cases). This shows that a scheme based on the 
GATE data would overestimate rain rate in the nidlatitudes, 
and additionally, if one uses a constant adjustment, there 


would be a large variance in ths rain anounts. Table XI 
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shows the comparison for GATE and Montreal rates of precipi- 
taticn. Wylie (1979) found that the precipitable water par- 
alleled the difference in the average rain rate with some 
minor excepticns when the troposphere was exceptionaily sta- 
ble. He decided to use Simpson and Wigqert's (1969) one-di- 


mensional cloud model to assess the effects of stability and 


to possibly improve his results in the mil-latitudes. 


Roinrote (mm h^!) 





280 270 260 250 240 


Infrared contour level (K) 


Figure 30. Rain Rates Averaged Over Cloud as Measured on 
IR Images for Various Cloud Threshol 
Temperatures (from Wylie, 1979) 
This model simulates cloud growth by using an 


ascending 2 km bubble. The modal incorporates cloud physics 


processes in the development of precipitation (Simpson and 
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TABLE XT 
Estimated and Observed Precipitation and Stabilitv Add 
4 o aT > ا‎ -— - ust = 
ment Factors (after Wylie, 1979) d 1 
(Montreal/GATE) Precipitation 
Precipi- 
Observed table Precipi- Model Adjusted 
rainrate? water table precipi- satellite Radar Error 
Case (mm h^!) (ram) water tation (10* m) Bias factor 
GATE 2-2 54 — =p Ba د‎ — — = 
Montreal 
29 June 1-0 38 0-7 0-5 24-6 26 5 9. l-1 
2 June 0-9 13 0-4 0-3 8-9 28 1-0 10 
26 Sept. 0-9 27 0.5 0-3 30-8 40-3 0-8 1-3 
l June 0-4 27 05 0-3 10-4 74 4 1.4 
20 Sept. 0-4 27 s 0-0 M2 14-6 0-2 4-6 
16 Sept. 0-1 22 0-4 0-0 7-0 7.2 1-0 1-0 
Average 14-1 17-5 0-9 1-7 
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NUJggert, 1969( Only one bubble ascent was used although 
normally the bubble would ascend many times. AS a result, 
the rainfall predicticn could not be used directly, but was 
nevertheless used as an adjustment to correct the S Tol 
Gama for stability. 

Table XI shows the results of Wylie's (1979) study: 
the corrected precipitation of the satellite closely matches 
the radar measured amount o£ rain. He noted that the over- 
all bias was small, and therefore by usin a stability cor- 
rection (calculated from the one-dimensional cloud model), 
the technique to measure precipitation by satellite in the 


tropics could be applied £o nid-latitude areas. 


r 
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III. SPADS CLOUD MIDEL DESCRIPTION 


A. INTRODUCTION 

The cloud analysis model adapted for the SPADS interac- 
tive computer is a combination of procedures drawn from sev- 
eral previously tested techniques. Tha mainstay cf the 


model is Liljas' (19812) parallelpipedon (Fig. 9 and 10) 


th 


classification. Added to this are a texture test for ċis- 
criminating some cloud types, a non-linear least squarss 
enye fit to discern cirrus and altostratus, Harris and Bar- 
Mees (1978) cloud amount calculation, part of Reynolds and 
Vonder Haar's (1977) cloud top temperature calculation, and 
Liljas' (19812) qualitative precipitation intensity model. 
The following sections describe these parts of the SPADS 
model and explain why they were selected. The last section 


of the chaptez summarizes with the flow chart of the model's 


computer progran. 


Eee CLOUD TYPING 
The apprcach selected for cloud typing was developed and 
tested by Liljas (1931a). Lt is a @ültispectrcal nethoG 


that uses three channels to type clouds and differentiate 
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between land and water. The SPADS cloud model uses two 
Channels, thereby becoming bispectral (the near-IR channel 
was not used due to non-availability in GOES VISSR data). 
Fig. 10 depicts the bispectral model iacorporated in the 
SPADS Cloud Model. Both the visual and infrared thresholds 
Eun in Fig. 10 were converted from Tiros-N Advanced Very 
High Resolution Radiometer (AVHRR) to GOES VISSR data 
counts. Liljas! visual data brightness counts ranged from 0 
to 255 whiie the SPADS counts range from I to 63. A four to 
one linear reiationship exists in the VIS channel and there- 
fore a direct conversion cf thresholds was performed. The 
IR thresholäs were converted to temperatıres (K) using an 


August 1980 TIROS-N AVHRR temperature calibration tabi 


(b 


| ۰ 


This gave a rain cloud threshoid of about 251 K which is in 
agreement with the Griffith/Woodley (Griffith and Woodley, 
1981) anā the Liljas (1981a) train cloud thresholds. These 
temperature thresholds were converted to GOES IR brightness 


counts. The transformed thresholds used in the SPADS cloud 


medel are depicted in Fig. 31. 
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One of the reasons for selecting the Liljas model 25 
that it was tested or an interactive computer similiar to 
SPADS. In ao dase lon, the test area selected by Liljas 
(1981a) included land and ocean areas along with many 56 
of clouds associated with mesoscale and synoptic scale fea- 
tures. À complete description with color pictures is in 
Liljas'( 1981a) technical report. rada tion, a combined 
approach using texture and bispectral counts was suggested 
by Liijas (1981a), Harris and Barrett (1978) and Reynolds 
and Vonder Haar (1977). The U.S. Air Force's 3DNEPH program 
used both these approaches, but documentation of thresholds 
and tests was not available to us. Thus, the 3DNEPH progran 
could only be used as a qualitative guide. 

Mee GOES IR resolution is approximately (Un mi. at the 
subsatellite point, but fcr the visual channel a higher res- 
Bon of 0.5 n mi. is available. This VIS resolution 
affords an cpportunity to get a better astimate of cloud 
cover. All the studies reviewed highly Tecommended che use 


an array of brightress counts to give an improved repre- 


th 


٥ 
Bean or the cloud pictures. The 0.5 rn mí. resolution 
allowed an 8X8 array of visual pixels par IR pixel. For 
each gri space an average visual brightness count is 


mmenjated (Eq. 1, Tabie II). 
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At the same time as the visual average brightness is 
being calculated, the standard deviation of the cloudy vis- 


ual counts per grid srace is deternined. The standard devi- 


ation (Eq. 4, Table I) is the statistic used to represent 
texture. As suggested and tested by Harris and Barrett 
(1978), the standard deviation is used to discriminate 
between stratiform and cumulus humulus, small and large 
cumulus congestus, and altostratus and cirrus. ThesinzeiaL 


standari deviation threshold values for this discriminant 
analysis were approximated by using information in Fig. 6 
and 8. An additional test was used first for “he cirrus/al- 
tostratus decision. A second order polynomial least squares 


fit was applied to the cirrus (c) points and also to the 


altostratus (d) points in Fig. 10. The calculated curves 
are depicted in Fig. 31, with tne associated equations in 
Table XII. Entering the equations with the average visual 


brightness, tne resulting IR value is compared *o the actual 
IR to determine cloud type. If the calculated value falls 
between the <wo curves then an additional test, using the 
Scandazi deviation tareshold, is used to discriminate cloud 


type. 
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(e, CLOUD TOP TEMPERATURES 

The Reynolds and Vonder Haar (1977) - study is used as a 
basis for estimating cloud top temperatures. TROIL 
calculation is based cn two equations (EJ. 1 and 2, Table 
IV), one for cloud amount and the other for the cloud top 
temperature. The resolution of the data used was U n ni. 
but th authors’ analysis suggested ar improved resolution 
would give better results. As nentioned above, this led to 
ene choice of 0.5 n mi, visual data in this study. 

The first step in determining the cloud top temperature 
is to calculate the amun* of cloud per IR pixel. The cloud 
decision was based on comparing the visual digital count of 


ee in the 8X8 grid <5 ano cloui +hreshoid. This 


(LU 


threshold value (initially 5 VIS counts) was derived fron 


Fig. 31 (after Liljas,  1981a) where it is the lowest VIS 


count to have clouds. A clcud threshold test was also used 
00 Harris and Barrett (1978). Their threshold value was 


derived by comparing the visual image with a computer pro- 


ensity slice, but neither the value nor the formula 


p 


duced 


Q 


Meor ain the value were included in their report. The U.S. 
Peereerorce’s 3DNEPH program based the cloud decision on the 


amount of visual grayshade difference between a 3 n ie 
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TABLED CET 


SPADS Çloud #odel's eye 1 Nec: Least squares 


pe 
Fquations, and Stan ard De? lation Valu ) 


Gi 


Cirrus (1)* - Cs 

a. y > 158.094 - 6.6076 x 4 4.159 x^ ** 
b. g < SIG3 

Altostratus (2) - As 

Bey 215,24 7.5240 x + .0395 x kk 
boo » STGJ 


* The number in the parens is the computer cloud type identifier. 


co A IR digital value 
x = average Visual digital value 


Fog/stratus (3) - ST 


Thick Fog/stratocumulus (4) - Sc 
ENS SIG 1 


Cumulus humilis (5) - CuHu 


O > SIG 1 


Small cumulus congestus (6) - CuCong 


O < SIG 2 


Large cumulus congestus (7) - CuCong 
g > SIG 2 
Nimbostratus/multilavered (8) - Ns 


Cumulonimbus (9) - Cb 
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pixel and a background brichtness (Fig. 15, Tefer to Chap- 
ter II for more complete details). This program uses an 
extensive glcbal backcround brightness file that is updated 
regularly (Fye, 1978). The development and maintenance of 
such a file for SPADS is not feasible at this time. 

The average cloud amount is determined by summing the 
number of cloud decisions in the grid space and dividing by 
Pees cOtalenumcoer of pixels per grid space. This calculation 
follcwed the Harris and Barrett (1978) Eq. 3 (Table 1). 
MARS cloud amount fraction is then used in the Reynolds and 
Wonder Haar (1977) equation for cloud top IR, radiance { Bq. 
D, Table IV). One of the reasons this cloud amount 
eal culation is used instead cf Reynolds and  Vonder Haar's 
emat on was that, in the latter, three aiditional constants 
would have to be calculated and verified, which is beyond 
Me Scope of this thesis. The authors felt that it was 
essential in using their cloud amount equation that these 
constants be well determined (Reynolds and Vonder Haar, 
1977). 

Mme razameter of the cloud top radiance equation to 
be determined is the clear area IÈ spectral radiance. Fleet 


Numerical Oceanography Center (FNOC) analyses of surface 
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temperatures  (K) are obtained and a representative 
temperature profile matching the location of the IR image 
center point is selected. These temperatures are converted 
to IR radiances via Planck's function and used in the cloud 
top temperature equation. 

Emissivity of the clouds also must be determined. 
Reynolds an Vonder Haar (1977), among others, found that 
the emissivity is nearly constant (approximately equal to 
Moye cor all clouds except cirrus, for Which it is approxi- 
mately 0.55. In the program, cirrus clouds are identified by 
the cloud typing and standard deviaticn nodules so that the 
lower emissivity can be used in the cloud top temperature 
equation. IT-T anvenssssvtv Of 1.0 was used for 
cumulonimbus and nimkostratus based on the assumption that 


these clouds do net allow any surface radiance through. 


The cioud top radiance (Ncld) is calculated using 
Reynolds and Vorder Haar’'s (1977) equation (Eq. 2, Table 
BTI) TST CEN 6025385 52562 ths  surtace “emperature field 


obtained from FNOC (cenverted to an IR radiance via Planck's 


0 


Bu, the emissivitv (citrus or non-cirrus), the satellit 
Memradtzance (Ml) value, and the £rac-tional amount of cloud 


uu). Once Ncld is calculated, a conversicn to a cloud top 
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temperature is made via the inverse of  Planck's function. 
This temperature can then be compared to the latest FNOC 
upper air data analyses to find the corresponding height in 
millibars (mb) of the cloud top. The upper air data are 
available at *en standard levals consisting of 1000, 900, 


oo, 700, 500, 400, 300, 250, 150, and 199 mb. 


D. ETUC ar TATION INTENSITY 

mnie precipitatior intensity identification techniques 
selected is the Liljas' (1981a) method. TE 1S an ev tens O 
of his cloud ‘typing sze*thod (see Chapter II) and therefore 
easily fits into the overall SPADS cloud model without add- 
ing much time to the process. dylie's (1979) cloud model, as 
ENcussed 2» Chapter II, used the Grifz£ith/Woodlev (Griz£zfi-*n 
IEEE. 1973) technique with a stability correction which 


Bee good indication ct rain rates, OR لاه اه دمه‎ gE- k ENS 


x 
e 


that the G Zith/Woodley technique is based on the develop- 
ment and identification cf convective cloud systems! life 
history through use cf radar echoes not available to this 
study. Another problem wich using this scheme is the 
computer run time it would add to the SPADS cloud model. 
BEEhough time 2S an important factor, it is not the only 


reason for selecting the Liljas method. Liljas  (1981a) 


reported good success with this technique, Obetelmiag dood 
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indications of relative precipitation over land and sea 
areas as verified by synoptic reports. He felt that this 
technique would be useful for sea areas where no synoptic 
observations were available. In addition, relative precipi- 
tation from mesoscale cloud clusters was detected between 


© مس ه 
a‏ 


ENUSUtic observations. This identificatio: 


b 
Ui 
N 
ny 
0 
tà 
(D 
tl 
للم‎ 
U) 
23 
O 
fu 
} 4 
i 


KS Eo. use in the SEADS cloud model by reducing the aunber 
BEDreci.pitation intensity catsjorizs fron six (see Fig. 9) 
Bes. Fig. 32 depicts the new Categories of precipita- 
tion types representing light, moderate, and heavy rainfall 


and Table XIII includes the equation for the lines s 


(U 


parat- 
ing the categories. whe light rainfall category includes 


the requirement that cloud top temperatures be less than 


Bo! K. 


Es LLOJ CHART 
The SPADS Cloud Mcdel (Fig. Hun "Degıns wich obtain!n 


ELM n--2:al ¿nirared and visual data fislis from the GOES 


tii 


sateilite and conventional surface and upper air data from 
IOC. The IR and VIS satellite data are collocated both in 
space and tine. The upper air and suriase temperatures ara 
E-Nected to correspond =o tne I,J position closest to tne 


center point of the image from the latest FNOC analysis on 
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the same date as the satellite pictures. E, ES ۹5 
cloud model computer program (see Appendix A) is run. 

As shown in Fig. 33, the first step of the program is to 
calcuiate three values based on the 8X8 grid of VIS pixels; 
the average VIS brightness, Standard deviation, and cloud 
amount as discussed in Section 2 of this chapter. Once the 
values are obtained, the cloud typing module (based on 
Liljas, 1981a) is called. Two tests are performed within 
this mciule; one is based on comparison of the IR and aver- 
age VIS counts and the other supplements the first by using 


texture (standard deviation) to complete the cloud typing 


وا 


task. If nimbostratus or cumulonimbus clouds are ident 


pcd, the precipitation module is calied and it 


fu 


etermines 


the intensity (qualitatively) of the precipitation. 


tJ. 
00 


The cloud top temperature module, in Watch, Ce 
temperature and the height of the clouds are calculated, is 
pcd for all cloud cases. The FNOC surface temperature 
value is used in the cloud top temperature calculation, 


while the  FNOC upper air data are sd in the height 


Es 
{A 
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esleulation. 
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TABLE XIII 
SPADS Cloud Model Precipitation Intensity Categories 


Equation suml = IR radiance + average VIS brightness 

No rain (d) suml < 184 

Light rain (1) 184 « suml « 195 

Moderate rain (2) 195 s sumPbos 224 

Heavy rain (3) 224 < suml 

The final step is processing the results. Cloud types 


ELE U- XII), —prec2pictacion inters2ta2es (Table XIII), cloud 
IE 'eemperetures (X), and cloud top heights (mb) are stcred 
in a permanent file fcr future display. The average visual 


Meegncmess, Standard deviazion, and the amount of cloud that 


or reference and 


Hi 


correspond to each IR pixel are printed 


validation purposes. 


on 





INPUT 
GOES IR Satellite Image 
GOES VIS Satellite Image 


FNOC Temp (SFC) 
FNOC Temp (UA mandatory) 














CALCULATE (8x8 sq. VIS) 
l. Average VIS Brightness 


2. Standard Deviation 
3. Cloud Amount 


DETERMINE CLOUD TYPE 
TR % AVG VIS -Thresholds 
. TEXTURE (standard deviation) 















mispeetral Calculation 
of Cloud Top, Temp 
and Height 







Determine Precip 
Intensity 
(Ns, Cb) 


Eure 33. The SPADS Cioud Model Generalized Flow Chart. 
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IV. THE MODEL TEST 


A. INTRODUCTION 


A test was conducted to datermine if the SPADS Cloud 
Model concepts were valid. The test was designed: (1) to 
capture a small area cf coincident IR and VIS satellite data 
in the mid-iatitudes and to obtain FNOC fields fcr the same 
date-time group and location as the satellite data, (2) to 
run the SPADS Cloud Mcdel using these data, and (3) to ana- 
lyze the results of the model to determin2 the validity of 
the modei's analyses. The following sections describe +he 
input *o this test, test procedures, results, and subsequent 
iterations and adjustments made to improve the model's 


performance. 


Bo CATA 

GOES IR and VIS data wera captured and archived on 
November 9, 1982 at 1315 local time, for a site centered ax 
location latitude 44°N, longitude 1410W. The corresponding 
FNOC data fields were also saved for the same date from 1200 
GMT data. Fig. 34 is the GOES IR image from which an area of 
apprcximately 460 X 8460 n ni. was selected a analysis. 


MWA-Warrorded 4)9G IR pixels for evaluation. Fig. 35 is the 
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IR satellite picture of the selected arsa with the center 
point matching that in Fig. 34. The IR data resolution was 
approximately 7 n mi. at this latitude and 4 n mi. at the 
subsateilite point. The visual data resolution used was 0.9 
nome. 10-5 Enz, at subsatellite point). The visual area 
(see Fig. 36) is the same area as the infrared giving an 8X3 
visual grid of pixels (64 points) per IR pixel. 


The FNOC vertical temperature profil was obtained -o 


(U 


match the center point cf the  Á satellites pictures and the 
corrsspanding heights (mb) are listed in Table XIV. Only one 
set of temperatures was used because the resolution or the 
FNOC fields (63 X 63 field with 381 kn mesh length) is 
approximately the same aS the area of the satelii-e 


pictures used in this test. 


C. mest DROCEDURES 
Men ollow ng list Is the set of procedures for this 
test of the SPADS Cloud Model. The steps taken were =o: 
eC Capture coincident VIS and IR satellite images. 


2. Receive surface and upper air temperatures from 
EOC. 


3. Manually analyze images Tor cloud boundaries and 
types. 


e Select the temperature profile closest to sataliite 
image 


e ct 


i 
a 


Q (P 
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5. Run SPADS Cloud Model conputer program. 
6. Output results in a format for analysis. 
7. Compare manual tc computer output. 

8. Adjust computer model as appropriate. 

9. Repeat steps 5, 5, and 7. 

The manual nephanaiysis was completed by Dr. C.H. Wash, 
Assistant Professor in Meteorology at the Navas Postgraduate 
School. He used both the IR and Visual images (Fig. 34, 35, 
and 36) to discern the cloud boundaries and types seen in 
Fig. 37 (see Tabie XV for symbol definition) 

ME Ag the first computer Eun, a cloud 5525526810052 5 
usual digital. counts was used. This valua corresporded to 
the lowest visual -hreshold that could hav2 clouds (see Fig. 
33). After analysis cf the first run, however, a threshold 


of 20 counts was selected. This adjustment accounts for our 


Ha 


Certainty in the GCES VISSR visible channel calibrazion, 


as compared <o that for the TIRƏS AVHRR visible channel. In 


in 


normal use, VISSR is a relatively qualitative tadiometer, 


and the engineering data needed +o convert values to high- 


ry 


quality calibrated tadiances were not readily availabie. 
Two sets of standard deviations were used and are presented 


Enable XVI. 
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TABLE XIV 


OC Temperature Profi ile or 1290 GM Nove poan 9, 1982 at 
E EH lat Den de 458 Ny ion m qd 14 19 


Standard level (mb) Temperature (°K) 


100 A 
150 214. 
200 216.8 
250 221.8 
300 O 229.] 
400 244.6 
500 - 254.3 
700 267.7 
850 275.5 

1000 283.4 

1010 (SFC) 289.7 
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NS 

Multiple Mid-Cloud 

Layers AS 
Possible Ac 





٩ 37.  Menual Clcud Analysis cf Fig. 35 and 36 (see 
Table XV for the symbci definitions). 
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چ سے 


Nimbostratus - Multiple layers 


Altcstratus 
Altocumulus 
Cirrus 
Cumulus 
Stratus 
Stratocumul us 


Broken 


ipai 


SIG (i) 
SIG 1 
SIG 2 


SIG 3 


un RESULTS 


TABLE XV 


Definitions of Symbols Used in Pig. 37 


NS -ML 
AS 
AC 
Ci 
Cu 
SE 
Se 


00 


TABLE XVI 


Standard Deviation (SIGi) Values Used for 
KIT EOI hoJ they 


T 1 and 2 (see 


are used 


Test 1 Test 2 
4 2779 
20 20 
3 5 


Mae Following subsections discuss the results of the two 


computer runs as 


compared دخ‎ 


the manual analysis. The 


results are discussed by the type of product and run number; 


101 





where the  prcducts are  cioud type, Przeıpıtation, cloud 


— 


amount, cloud top temperature and height, and texture 
(standard deviation). 
1. Run 1 
a. Cloud Typing 

The overall general patterns of the cloud types 
(see Fig. 38) compared well to the manual analysis (see rig. 
B7). The observed cloud structure is most easily discussed 
by considering the images as four quadrants; northeast (NE), 
northwest (NW), . southwest (S4), and southeast (SE) with 
north being the Top of the image. 


mid-level clouds such as 


in 


The NE quadrant ha 
EBMEOStTratus (see Fiq, 37). Tha  SPADS cloud model azd well 
here indicating predominatly mid-level clouds (see Fig. 38). 
The model also matched the manual analysis in the NW quad- 
mant. This area has high cloud tops of nimbostratus and 
embedded cumuloninbus. 

The SW quadrant has a mixture of mid-level and 
multi-layered thick clouds and again the cloud model is in 
agreement with the manual analysis (see Fig. 37 and 38). The 
objective analysis shows mora detail than tne manual 


enalysis because the computer program is working pixel by 
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pixel while the manual analysis can only resolve groups of 
pixel defining cloud areas. 

In the SE quadrant, the SPADS cloud model could 
not resolve the difference between low-level (stratus) and 
mid-level clouds. The boundary between these is located on 
the upper border of the SE quadrant. Another discrimination 
problem is between stratus and cumulus in the lower right 
quadrant. The amount of cumuius humulus that seems apparent 
in the visual satellite picture (ses Fig. 36) and the manual 
analysis (see Fig. 37) did not appear in the computer 
analysis (see Fig. 38). This problem is related to the tex- 
ture threshold tests, as bozh stratus 2534 cumulus humuius 


have the same cloud type thresholds. 


in 


EN Precipitation 

Biesbatraernszarsrannenwensity (see Fig. 39) in 
the NW and SW quadrants are based on the premise that the 
Breker and colder, tnerefore orighter, clouds (compare Fiq. 
EO Fig. 35 and 36) have a greater probability of rain and 
DE SUuer in*ersity of rain. Along with this, is the knowl- 
edge that the approximate threshold for rain clouds is 
colder than -22 °C as ۳۹31 : Barrett'and Martin, 


1981). Both computer runs gave the sane agreement with the 
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Figure 38. Cloud Types as Analyzed by ths SPADS Cloud Medel 
in Run I (see Table XII for numeric 
definitions). 





nanual analysis because the areas of rain clouds did no- 

change. At this time, no surface data are available for 

ground-truth verification of the rain patterns. 
Co Cloud Amount 

The cloud amount for the first run was 100% for 


١ 


the entire analysis area. The visual, IR, and manual pic- 


LY 


tures ail shcwed clear and scattered cloud areas in the § 
quadrant. In the manual analysis (see Fig. 37), this quad- 
rant has these sections marked Cu, St/Sc,. and 2 AS (see 
Table XV for symbol definitions}. Further research into the 
cloud threshold for GOES satellite imagery, combined with 
analysis of the high resolution visual data used in this 
test, suggested a better cloud +hreshold value would be 20 
visual counts (Muench and Keegan, 1979, and see discussion 
Me ection C). The results from using this thrsshold are 
discussed in tne Run 2 section. 
d. Cloud Top Temperatures and Heights 


The cloud top temperatures )5 2 40) and 


(b 
ل۱‎ 
I^ 

Q 
e 


heights (see Fig. 41) look reasonable for the cloud types 
BEmtULried. The NE quadrant had temperatures ranging from 
260 to 250 X and a telght or approixmately 500 mb, these 


matching the type of mid-level cloud. The NW quadrant 
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Eure 39. 117 ce ve Proci pi at Lor In 
of the SPADS Cloud Model (s 
pume: e definitions. 
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clouds are mostly nimtostratus with high-level, cold cloud 
tops. The height is approximately 300 mb and tne 
temperature is 230 K. Again, agreement is good. The mixed 
clouds in the SW quadrant also have mixed height (mid-level 
and high-level) and temperatures which compare well (see 


Fig. 40 and 41) 


tad 


E SA CSS 15ص‎ Ban nom. tom gol so 
700 mb, define the area reasonably well (see Fig. 41). The 
division for the low and mid-level clouds (see Fig. 37) 


appears as the 700 mb contour line (see Fig. 40). The cloud 


temperature field (see Fig. 40) defines the division of low 


and miá-levei clouds less clearly, but partly cloudy areas 
|| Kk contour)aze indicated. BNEWEJEENS and altostzarus 
SONOS are too cold and high, aS compazsa to the  FNOC 


temperature profile. This error appears to be in the appli- 
cation of the bispectral radiance equation (Eq. 4, Table IV) 


from Reynolds and Vender Haar (1979) in partly cloudy 


regicns. The ccefficients in that equation depend strorgly 
on the radiance calibration of the radiometer used. The 
assumed calibration ccefficsints used in this study are 


Cleariy inadequate and must be improved through empirical 
udation of GOES VISSR calibration asing ground truth 


data. 
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Cloud Top Heights (mb) from Run 1 of the SPADS 


Cloud Model. 


Figure 41. 





e.  Clcud Texture (Standard Deviation). 

The clo ud texture (standard deviation) 
calculation results are presented in Fig. 42. Overall, they 
show an expected pattern; mixed and broken cloud areas (in 
the SW and SE quadrants) have the highest variation and uni- 
form stratus type cloud cover areas in the NW and NE quad- 
rants) have the smallest variation (compare Fig. 36 and 42). 
The highest deviations (approximately 4.0) were located in 
the SE guadrant in the areas discussed previously under 
cloud amount. After contouring, there appeared to be a defi- 
nite boundary between the highly variable region (SE quad- 
rant) and the nearly uniform areas (NE quadrant, see Fig. 
42). The SE quadrant had been manually typed as cumulus, 
with broken areas, and stratus/stratocumulus clouds (see 
EET 37). After evaluation, the standard deviation threshold 
value for discerning cumulus humulus was adjusted along with 
the value for altostratus. There did not appear to be obvi- 
ous patterns of different size cumulus congestus and there- 
fore this standard deviation threshold was net adjusted. 
The standard deviation definitely shows promise in 


discriminating between cumuloform and stratus cloud types. 
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The results from *he second run are presented below 
with the revised cloud and standard deviation thresholds. 
Only regions where changes from Run 1 occurred are discussed 
in detail. 

a. Cloud Typing 

The cloud typing results were changed somewhat 

because of the new cloud threshold and standard deviation 
thresholds (see Fig. 43) especially in the lower right quad- 
mant. More cumulus humulus and altostatus types were deli- 
neated and the patterns appear closer to those in the manual 
analysis. 

DO Ciud Àmount 

The cloud amount patterns (see Fig. 44) are nore 

in line with the manual analysis and what can clearly pe 
seen on the satellite images. It appears thatthe cloud 
enreehnld 0£ 209 is better, but this may require additional 
study. 

oe O Top Temperecare and Hezght 
Tie BOBBLE in obzainıng the ccrrect cloud top 


temperature and height in a partly cloudy situation using a 


U 
ct 
(D 


LU 
F^ 
00 


di 


0 


bispectrai technique still remained. This problem 


Uu 


mom uncertainty in calibration of the GOES VISSR. Ihe VISSR 
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infrared channel is not actively calibrated and extensive 
ground truth is required to accurately interpret its data in 
terms of temperatures (Maul, 1981). In the present program, 
Pee erature calibrati كك‎ would anter into the coefficients of 
the cloud height Eq. 4 (Table IV) and directly into the 


cloud top temperatures. Time was not available in this study 


1 


Te 


سا 


to make this adjustme 
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Ae DISCUSSION 

This chapter summarizes the thesis accomplishments, the 
problems encountered, and suggestions for future research. 

1. Accomplisnments 

The first accomplishment was to complete an overview 
of cloud and precipitation studies with an emphasis given to 
those using an interactive computer system. Next, useful 
techniques and ideas were selacted from these studies and 
used to design tke SPADS cloud and precipitation model. This 
SPADS Cloud Model was then coded and established on the 
SPADS systen. Finally, two tast runs on a test case were 
completed with good success. 

NIN dN INOUE, detailed intormation was produced by 
the SPADS Cloud Model through concurrent analysis of high 
mecolution visual and infrared digital satellite data. Over 
a maritime region, the SPADS Cloud Model estimated various 
mould -ypes, precipitation intensity, cloud amounts, and 


cloud top temperatures and heigkts. 
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roblems 


N) 
è 
0 


Several problems were encountered. The cloud amount 
threshold needed adjustment and two tests of different 
thresholds were acccmpiished. The tareshold of 20 was 
satisfactory but furtter verification is suggested. 


y 


Another zrea cf concern is the threshold values for 
the texture tests. Two thresholds were adjusted in the tests 
but all three need further refinement. Ihe separation test 
(the least square fit equations) co rMeSrrus and altostra.us 
could not be substantiated in this data sez and will aiso 
need verification in future runs. 

An additional problen Pen the cloud top 
temperature and heign+ calculations. ise Sscusscq SEO; 
the temperatures are too cold and the heights too high in 
eertain partly cloudy areas. The bispectral equation to 
obtain these values needs further study. 

Ali tnree of these problems are related, at least in 
part, to the problem of inter-calibration of AVHRR and GOES 


ISSR visible and infrared da-s. 
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The following are recommendations and suggestions 
for additional research. The SPADS Cloud Model should be: 


a. converted to use temperatures for infrared thresholds 
and albedoes for visual thresnoids to allow easier con- 
version between different satellites, (polar and geosta- 
tionary) and seasons (this will require development and 
addition of sensor calibration tables and modules to 
the SPADS Cloud Model). 


converted tc a variable resclution so that larger areas 
could be analyzed in the same time frame but also keep 
the high resolution capability. 


c. tested and verified in various regions by using surface 
observations, _ Tawinsondes, 31272 puattorm 
observetlons, ship observations, radar,  and,if avaiia- 
ple, the U.S. Air Pcrce 3DNEPH cloud analysis. 


Digital satellite data are currently underutilized, 
except perhaps for sea surface temperature and temperature 
sounding applications. Most visual and infrared satellite 
data are treated as images and are stili analyzed manualiy. 
Now, with the development of SPADS and similiar systems, new 
satellite products can be produced in regional forecast 
Sites based or digital satellite data. This thesis illus- 
trates the detailed analysis of the digital satellite infor- 
mation that can be received by the user ina short (l-2 
hours) time frame. Specifically, this work was directed to 
using the SPADS comcuter system toward process satellite 
visual and infrared digital data and to produce cloud and 


Precipitation inZornation. As SPADS is placed in the field 
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reat sites, this infecogation will, for the first tine, 
be made available in an operational setting instead of being 
available only to meteorologists at the major satellite cen- 
ters. In the future, this 2nd other applications wiil sub- 


stantially aic the operational forecaster. 
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APPENDIX A 


SPADS CLOUD MODEL COMPUTER PROGRAM 


CLOD -- IS DRIVER/MAIN PROGRAM 


BIS ZEROGRAOMZANAETZES VIS AND IF TOGETHER TO DERIVE TYPES OF 
ESOUDSSCDOUDOSHETGHTSS.CLOUD AMOUNT, CLOUD TOP TEMPERATURES, 
BNIDMRELATIVE INTENSITY OF RAIN. 

PROGRAMMED BY LT. C. NELSON 

INITIALIZATION 


COMMON/A/ML(128) CIRRUS,K,KI,ISEL 

COMMON/B3/CTVXPE(64) PRCIPI(64) ,MSAVG(64),1 
COMMCN/C/SIGMA(64) ,5 16151162 , 3 

GO نين‎ 15 1007 )64 ,CTEMP(64) ACDAVG(64),E1,E2,ALAT,ALON, Jl, 


BBD), ES 


COLMON/F/ZP.,XCEN, XCENSELIPJAL, 3ASLOHĦ,LVL(11) 
COMMON/G/STDP(11) 
E PNIS? ONDUSMSAVGISIGMASEL,E2,E,ACLD, ACDAVG, AMCLD, SUM, SUMSIG, 


*HEIGHT,N,CTEMP,SIGl,SIG2,SIG3 


INTEGER CTVPE,CIRRUS PRCIPL,M,0,P,L,PCLDI,PCLD2,0,ISEL,MS(256,3), 


*VISPN(7) , IREN(7) ,MS1(8) 


Il=1 
ESEL-0 

M20 

N=0.9 
E120.9 
E220.55 
ESSO. 
SUM=0.0 
SUMSIG20.0 
AMCLD=0.0 
CIRRUS=0 
ACLD=N.0 
P264 
SIGI=2.5 
SIG2=20.0 
SIG3-21.5 
PCLD1220.0 
02512 
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170 


1:50 


i— we) 
A r3 
سم‎ TO SN 


190 


MENU E 


TYPE "4155" 

TYPE "THIS PROGRAM IDENTIFIES CLOUD TYPES" 
TYPE "ADDITIONAL CLOUD DATA CAN ALSO BE OBTAINED" 
TYPE "THESE ARE - " 

TYPE "0 NO ADDITIONAL INFORMATION" 

TYPE "] PRECIPITATION INTENSITY" 

uu "ONGLOUD TOP TEMPESATBRE" 

TYPE "3 CLOUD TOP HEIGHT" 

TYPE "4 ALL OF THE ABQUE" 

TYDE "8 EXIT PROGRAM" 

ACCEPT "ENTER REQUEST= ",ISFL 
EISEL. EQ 38)60 TO 80 
MEGESEIWEN.O.0OR.ISEL-E0:1)GO TO 131 

GNLL OPEN(A," TLVN.O" , 1, IER) 

IF(IER.EQ.1)GO TO 170 

TYPE "TLV FILE OPEN ERROR, RC: ",IER 

CALL RDBLK(4,1,!S,L, IZR) 

IF(IER.EO,1)GO TO 180 


TYPE "TLV FILE PEAD ERROR, 2C= ",IER 
CO TO 70 

00 1 I=1,11 
TLV(I)=FLOAT(US(1))/100.0 

CONTINUE 


gwen s( 12,626) (TLV(L),L=1,11) 
FORMAT(2X,11F5.2) 

FORMAT( S13) 

EPE "<15>" 

ACCEPT "ENTER IR FILENAME" 
READ(I1,2) IREN(CI) 

SEL OPEN(2,IREN,1,IER) 
MES SEO. 1)60 TO 190 

FR FILE =~ OPEN ERROR, RC= ",IER‏ قر ل 
meee. “TRY ACAIN?” IX‏ 

XU EO.1)G6G0 TO 181 
MT. EN -0)GO TO 70 

ACCEPT "ENTER VIS FILENAME" 
RESDCLI.2)VISFN(I) 


123 
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CALL OPEN(1,VISEN,l,IER) 
IF(IER.EO.1)GO TO 200 

TYPE "VIS FILE - OPFN ERROR, RC» ",IER 
NECEPT "TRY AGAIN?" ,IX 

IF(IX.EO.1)GO TO 190 

GO TO 69 


OPEN FILES FOR STORAGE OF CALCULATED VALUES 


CALL PEN(S, CALCID.0",3, IER) 
IF(IER.EO.1)G60 TO 201 

TYPE "DATA 1 FILE OPFI ERROR, AC= "LIER 
GO TO 69 

CNN. OPEN(7 , CALC2D. 
IF(IER.E0.1)GO TO 20 
TYPE "DATA 2 FILE OPEN ERROR, RC= ",IER 
GO TO 69 


DCO TÛ 
D 


DO LOOP TO CALL IN SERIES OF DATA 

DEDE S RUN PROMO TQ 511 

MNMEMEMARES AT 192,192 FOR IR TO MATCH VIS (COLLOCATED CENTER 
BOENT). IR I5 2 MILE E FIORE IED TO USE EVERY GIZFER VALUE 


I320 

12=0 

noo? TIs1. 126,2 
I3=(Il- 1)+191 
12=1+12 

DO 6 I=1,P 
HEIGHT(I)= 
PRCIPI(I) 
ACDAVG(I)s 
SIGMA(I)20.0 
CTVPE(I)20 
CTEMP(I)20,0 
CONTINUE 


1 


0.0 
ie) 
0.0 


DRIVER MODULE VHICH WILL PREPARE VISUAL DATA TO 35 USED IN OTHER 
MODULES. THE AVERAGE VISUAL BRIGHTNESS, THE STANDARD DEVIATION OF 
VISUAL 3RIGHTNESS, AND THE CLOUD AMONT IN 8X8 GRIN SETS ILL 3E 





C3 


IN 


Q 


204 


10 
P 


DETERMINES. 
IREDTGITALSVSEUE FILE 


CALL RDBLK(2,13,MS(1,1),1,IER) 
IF(IER.EQ.1)GO TO 203 

TYPE "IR FILE - READ ERROR, RC= ",IER 
GO TO 69 

CALL UPK(MS(96,1),MS(1,2),128) 

nn 100 I=1,123 

ML(I) =FLOAT(MS(I,2)) 

CONTINUE 

€120 

J1=(12-1)*3 

Pegi BOBEX(1,J1,1S,3,1ER) 
IF(IER.70.1)GO TO 204 

7YPE "VISUAL FILE READ ERROR, RC= ",IZR 
GO TO 69 

CONTINUE 

36 30 %251,128,2 


merece PROCESSES VISUAL DATA 3Y 9X3 GRID 


CONTINUE 

CONTINUF 

L1=(X1-1)*4+1 

90 11 0=1,3 

CALL UPKOIS(L1,0),'181,3 
50 10 L=1,2 
ECODPEDOSFLOAT(HSICL)) 
CONTINUE 

CONTINUE 

11-8 

20513 L-1,3 

DO 12 0-11 


[2 
A 


Fr 





kJ 


mm 


DOSS 


2 
13 


13 


30 


N=N+1.0 

SUM=SUM+MS2(L,Q) 

SUMSIGSSUMSIGT(MS2(L,0))**2.,0 

IF(MS2(L,0).GE.PCLD1)ACLD-1.0 
IF(4MS2(L,O).LT.PCLDL)ACLD=0.0 

AMCLDSACLD--AMCLD 

CONTINUE 

CONTINUE 
SIGMA(X1)5SORT((N*SUMSIG-SUM**2,0)/(N*(N-1.0))) 

MSAVG(K1)=SUM/N 

CALL CLOUD 

ABBAVG(EI)SAMGLDJM 
BEOOTYPE(K1).£E9.8, OB, CTYPE(KI) . EO. 2) 2221.0 
NSGNSEECEO.OQOR.ISEL.EO.1)GO TO 13 

CALL TEMHT 

LF(ISEL.EQ.S)CN TO 69 

CIRRUSE0 

N&0.0 

ACLN=N.O 

AMCLD=0.0 

SU: 120.0 

SUMSIG=0.0 

E3=0.0 

CONTINUE 


STORE ALL VALUES 
Zee STORED IN CALC] FILE =— CTYP",ePCIPI,CTEPP, 
HEIGHT(1-256 WD) 

emus IN CALC2 — HSAVG,SIGIA,ACDAVG OC 1-256 WD) 


n0 23 I=1,64 
MS(I,L)=CTYPE(T) 

J= 1+64 

MS(J,L)=PRCIPI(I) 

X5=J+64 
MS(K5,1)=IFIX(CTEMP(I)+.5) 

X4AzK 54-54 

MS(KA4 1)eI7IX(HEIGĦT(I)4.5) 

SI 2 = EE IMC MSAVG(T)+.5) 

BOSCO I OIFIXCIOO.O0fSIGMA(I)-T.5) 


126 





24 
25 
27 


E 
69 


70 


AS(K5,2)=IFIX(1N0.0%ACDAVG(I)+.5) 
MS(K4,2)=0 

CONTINUE 

K7=12-1 

CALL WRBLK(5,K7,MS(1,1),1,IER) 
Pu@uER. EO.) )GO TO 22 

TYPE "WRITE ERROR 5 RE= " IER 
GO TO 69 

CALL WRBLK(7 ,K7 p4S(1,2),1, IER) 
MBR SEO. 1)CO TO 26 

IE US ITE ERROR. 7 RE= "LIER 
GO TO 69 

EU 12 27) CHEIGHT(I),Is1,64) 
WRETDECI? ,25)( SIGMA(I),Is1,54) 
Bell? 27 )CCTEMP(T),T=1,6&) 
WRITE(12,24) (MSAVG(I),I=1,64) 
NRETECI2,24) CACDAVG(I),Is1,64) 
Emu 25) )(CTYPE(L),1=1,64) 
WRITE(12,25)(PRCIPI(I),Is1,64) 
Pompe 2.27) (MLCL) , 122,128 ,2) 
BOBMAT(/ 4(/,15(IX,F7.4))) 
FORMAT(/ ,6412) 

OMA TC), 16C1X,F7.2))) 
Dee. Pr CONTINUE?" ISTP 


.IF(ISTP.EQ.1)60 TO 69 


CONTINUE 
CALL CLOSE(l,IER) 
CALL CLOSE(2,IER) 
CALL CLOSE(4, LEX) 
CALL CLOS&(5, IER) 
CALL CLOSE(7 ,IER) 
STOP 

END 
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SUBROUTINE CLOUD 


ERRBEBETNEINGZODUER, ALSOZTHE PRECIPITATION MODULE TS CALLED 
FROM HERE WHEN PRECIPITATION CLOUDS ARE IDENTIFIED. 


COMMON/ A/ML(128) ,CIRRUS,K,K1,ISEL 
COMMON/B/CTYPE(64) ,PRCIPI(64) ,MSAVG(64) ,™ 

REAL ML,MSAYG, HEIGHT, CTEMP,SIG,SIG1,SIG2,SIG3 
INTEGER 1,CTYPE,PRCIPI,CIRRUS,ISEL 


Piscean LIAS CONVERTED GOES THRESHOLDS 


H2= 24.0 
H32100.0 
H4=113.0 
H52159.0 
H6=255.0 


THE FOLLOWING IF STATEMENTS CHECK 30TH IR (ML) AND VIS (HSAVG) 
VALUES AGAINST THRESHOLDS TO DETERMINE CLOUD TYPE. THE 
DEANDARDODEVIATION (SIG) OF THE 8X8 VISUAL GRID IS USED. 


TYPE ] AND 2 ---- CIRBRUS/ALTOSTRATUS 
ITE DIGTPOHASANDOMECTY).LE.HG.AND.MSAVG(K1).GT.5.0.AND. 
ASAS (1). .L2.23.0)60 TO 700 
TYPE 3 -=== FOG/STRATUS/CUMULUS 115 
IU NE TELA ا‎ HA SASDEMSAUGCKL).GT.13.0.AD. 
ENEAVGXI).LE.20.0)GO TO 900 
TYPE 4 AND 5 ---- THICK FOG/STRATOGUMULUS/CU'ULUS 5 
ال‎ CO.GT.HO. AMD. AL(R).LE.HA.AND.MSAVG(K1).GT.20.0.AND. 
EISSUGOKI).LE.36.0)GO TO. 600 
GO TO 601 
600 “=l 
GO TO 800 
TYPE 6 AND 7 ---- SMALL/LARGE CUMULUS CONGESTUS 
BIMIMMBEC(VE(K)GTAĠ AND, ML(E).LE,45 .,-AND,MSAVG(EI).GT.23.0.AND, 
XUSAVG(KI).LE.39,0)GO TO 605 
GO TO 606 
605 2 
GO TO 800 


r— 
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TYPE 8 -—-- NIMBOSTRATUS (MULTI-LAYERED) 
EUER GROH5,.ANDSMIEQOEOTESSHO. AND, MSAVGCKI) .GT,23.0. AND, 
*MSAVG(X1).LE.31.0)GO TO 910 
TYPE 9 —-- CUMULONIMBUS 
INE GTOHSCANDEMECROSEECHS AND. SAVGCK1).GT.31.0) 
*GO TO 915 
CTVPE(RI)20 
00 TO 8 
CALL STDDEV MODULE TO FURTHER DISCRIMINATE CLOUD TYPES 
700 M23 
300 CALL STDDEV 
GO TO 33 
LABEL TYPES 
900 M=4 
GO TO 800 
38 PRCIPI(K1)=0 
COTO 40 
910 CTYPE(X1)=8 
CO TO 970 
DES CTYPE(K] )=9 
BONO JIND 9 IVPES ARE CONSIDERED RAIM CLOUDS, SO GALL 
PRECIP MODULE. 
A SEL SEO. 0.0... SEL. EQ.<2. OR. ISEL.EQ.3)G0 TN 40 
CALL PRECIP 
40 RETURN 
IND 


pe 





MS OOO 


O 
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IC 


SUBROUTINE STDDEV 


STANDARD DEVIATION TEST MODULE HELPS DISCRIMINATE FURTHER THE 
MMC RISO LT SECIS CIRRUS FLAG FOR USE IN THE TEMHT 
MODULE. 


COMMON/ A/ML( 128) ,CIRRUS,K,K1,ISEL 
COMMON/B/CTYPE(A4) ,PRCIPI( 64) ,MSAVG(64) ,M 
COMMON / G/STGMACH4) ,SIGL ,SIG2 ,SIG3 

REAL *{L,MSAVG, HEIGHT,CTEMP ,SIG,Y,SIG1,SIG2,SIG3 
meer CTYPE, PRCIPI,CIRRUS, 1,1 


IS(M.EQ.4)G0 TO 813 
IF(M.EQ.L)GO TO 810 
IF(M.EO.2)CH TO 313 
CIRRUS VS ALTOSTRATUS, FIRST TEST : 
THIS EQUATION IS TOR LILJAS CONVERTED GOES THRESHOLDS 
Y=158.094-6 ..5076*MSAVG(X1) +4. 159%*MSAVG(K1)%*2.0 
IF(Y.LE.ML(X))GO TO 814 
SECOND TEST FOR CI/AS 
THIS EQUATION IS FOR LILJAS CONVERTED GOES THRFSHOLDS 
Y2115.24-0,5246* 4SAVG(1 )-0,0395s19AVG( 1) 0 
LF(Y.GE.ML(X))GN TO 316 
THIRD TEST FOR CI/AS 
EF(SIGMA(K1):GT.SIGC3)GO TO 314 
GO TO 815 
EMOMITE(SIGĦA(LIJELT.SIG1)6O 70 311 
TYPE 5 — CUMULUS HUMULUS 
CTYPE(X])=5 
0 10 320 
TYPE 4 — THICK FOG/STRATOCUMULUS 
311 CTYPE(X1)=4 
GO TO 820 
pis) LP (SIGMA(K1).LT.SIG2)GO TO 817 
CTVPE(RI)26 
TYPE 6 — SMALL CUMULUS CONGESTUS 
GO TO 820 
TYPE 1 -- CIRRUS 
314 CTYPE(K1)=1 
CIRRUS=1 





e 


e 


816 


817 


818 


D 
9 


GO TO 820 

TYPE 2 — ALTOSTRATUS 
CTYPE(KL)=2 

GO TO 820 

TYPE 7 -- CUMULUS CONGESTUS - LARGE 
CTYPE(K1)=7 

GO TO 820 

TELSTENALKT )SLT.SIGL)GO TO 819 
TYPE 5 — CUMULUS HUMULUS 
[DEC )-5 

GO TO 820 

TYPE 3 — FOG/STRATUS 
CTYPE(X1)23 

RETURN 

END 


131 











= 0 


١ 


toc mus CHAM 
con BEIM em I 

cot D. 

i DIN 

CPWNT سر‎ 


Lei LT 


wA 





OOOO 


ISO 


دی نت ده 


SUBROUTINE PRECIP 


THLS MODULE IDENTIFIES PRECIPITATION INTENSITY AREAS 
1=LIGHT, 2=MODERATE, 3=HEAVY, AND O-NO PAIN. 


COMMON/ A/ML(128) , CIRRUS, K,K1, ISEL 

COMMON/ B/CTYPE(64) ,PRCIPI(64),MSAVG( 64) ,M 
REAL ML,HEIGHT,CTEMP,MSAVG 

INTEGER X,CIRRUS,PRCIPI,CTYPE,SUM] 


ILIBJAS CONVERTED COES THRESHOLDS‏ ااال 


11=134 
12=195 
13=224 


ADD IR (ML) AND AVERAGE VISUAL (MSAVG), THIS SUM YILL BE TESTED 
FOR INTENSITY OF RAIN. 


SUMIEML(K)-FISAVG(KI) 

EHE UN CESDUSAND.SUMI.LE.I2)GO- TO. 43 
E GCITSIP2TAND.SUMI.LE.I3)GO TO 44 
TEE .CT,1L3)6C0 TO 45 

PRCIPI(KI)20 

GO TO 49 

PRCIPI(Kl)=1 

GO TO 49 

PRCIPI(KI)22 

GO TO 49 

PRCIPI(KI)23 

RETURN 


END 
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SUBROUTINE TEMHT 


THIS SUBROUTINE CALCULATES THE CLOUD TOP TEMPERATURE (CTEMP) 
USING THE IR RADIANCE FOR THE SURFACE (NCLR), THE AMOUNT OF CLOUD 
(ACDAVG), EMISSIVITY (E), AND THE MEASURED IR RADIANCE FROM THE 
SATELLITE (ML). CISMP IS THEN CORRELATED TO AN AREA 
REPRESENTATIVE UPPER AIR SOUNDING (UALV) TO DETERMINE THE CLOUD 
TOP HEIGHT (HEIGHT). 


CONBION/ A/MLCI28) , CIRRUS,K,K1, ISEL 
COMMON/D/'IEIGHT(64) ,CT! EMP( 6A) , ACWAVG( 64) ,E1,E2,ALAT,ALON,J1, 
*TLV(11),E3 
"a 5 
اا‎ HEIGHEISCTEMP,EL,E2,NCER,NCLD,ACDAVG,E,SIGI,SIC2,SIG3, 
*TLV,STIP 
INTEGER I,CTYPE,PRCIPI,CIRRUS 
PROASS iD? / 10100, 1000.07350,0,700.0,500.0,400.0,300.0,250.0,200.0, 
150.07 
E=0.0 
CHECK TO SEE IF THERE ARE ANY CLOUDS, IF NOT SKIP TO END 
1000000 0 00617 
CONVERT SURFACE TEMPERATURE (TSFC) TO NCLR 
meme (ll )).Gre242.0)NCLR=2.0*(331.0-TLV(1)) 
Be Tuch) .200242.9)1CL3=420.0-TLV(1) 
IS CIRRUS FLAG SET? IF YES, USE F2 FOR E. IF NOT, USE El FOR E. 
DEQCTISRUS.EO.1)CO TO- R50 
IS THE C3 OR NS CLOUD FLAG SET? IF YES USE E3 FOR El. 
fer (293.80.1.0)G0 TO 849 
E=E] 
GO TO 851 
849 E=E3 
CO TO 851 
850 EsE2 
COMPUTE CLOUD RADIANCE 
851 NCLOZML(K)-NCLR/(E*ACDAVG(K1) )--NCLR 
CONVERT TO CTEMP 
BESGICSDOEE.078.0) CTÉMP(K1)2331.0-'ICLD*0.5 
IF(NCLD.CT.178.0) CTEMP(K1)2420.0-NCLD 
WRITE(I2,853)NCLR,NCLD 
ECSSEDRMATC/IZ,5(1H*,E15.3) ,  NCLR, ICLD") 


D 
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Bl? 


313 


316 


314 
315 
60 


51 


IF(ISEL.EQ.2)G0 TO 60 


INTERPOLATE UPPZR AIR SOUNDING TO CONTINUOUS CURVE OF TEMP 
e U GERT, THEN COMPARE TO CTEMP TILL GET MATCH. CRECK 
MEONS TORS DOWN. TREN USE HATCH TO SELECT CORRECT HEIGHT FOR 
CLOUD TOP. 


DO 812 I822,11 

IAz13-I3 

Pere tow) toms tLV(11).GTSCTEMP(K1))GO TO 814 

NUR SEA .LPRETEMPEKI))GO TO 313 
IEDE,SOPCTEMPMI))CO TO S14 

CONTINUE 

IF(IA.EQ.2)1A23 
RHGSEXP(ALOG(STDP( IA-1)) -CALOG(STDP( L3-1) ) -ALOG( STDP( 1A) ) ) * 


BPEUGT FID-CTEMP(KI1))/(TLV(TAFI)-TLV(TA))) 


GO TO 815 

EX NU PERATURE IS ABOVE 100 MB" 
ARRDPIOUDONTENUE?'".TCONT 
IFCICONT.EO.N) ISEL=8 
RHG=STNP(IA) 

HEIGHT(K1 )=RHG+0.5 
RETURN 

GTEMP(£1)-TLV(1) 
CEEI SS TBC y 
RETURN 

END 





10. 


IEDSIEMAL DISTRIBUTION LIST 


NO. 


Defense Technical Information Center 
Cameron Station , 
Alexandria, VA 22314 


Library, Code 0142 
Naval Postqraduate School 
Monterey, CA 93940 


PDroresscr Acbert J. Renard, Code 6381 
Department of nS oco ay 

Naval Postgraduate Schoo 

Monterey, CA 93940 


PEOrEcssO> Christcpher N. K. Mooers, Zoda 68 
Department of Oceanography 

Naval a an, School 

Monterey, CA 93940 


Assistant Professor Cariylə H. Wash, Code 63Wy 
Department of Meteorolog 

Naval Postqradua*e Schoo 

Monterey, CA 93940 


Assistant Professor James L. Mueller, Code 684y 
Departmert of Oceanography 

Naval Postgraduate School 

Monterey, CA 93940 


Lt. Cynthia Ann Neison 
Department of Oceanography 
U.S. Naval Acaden 
Annapolis, MD 21402 


DSEG G C . 

Naval Oceanography Division 
Naval Observatory 

34th and Massachusetts Avenue NY 
nace naton, D.C. 20390 


Conmander 

Naval Oceanograpty Command 
ASE station 

Bay ot. LOUIS, MS 39522 


Commanding Officer 
Txal Oceanographic Office 
foe Station 

EY ST. Louis, MS 39522 


135 


Copies 





11 


22. 


i. 


EJ 


15. 


16 


17 


13. 


19. 


20. 


Conmanding 
Fleet une 
ee 


C 


fi 
SET 
A 93 


Commandirg officer 


Naval Ocean Research and Development 


Activit 


NSTL Station 


Bay St. LOUIS, 


Commanding Officer 
Naval Environmental Prediction Research 


BC 


Facili 
Monterey, 


A 93940 


CDR Don Hinsman 


Naval E mental Prediction Research 


Facil 
Monterey, 


“E 


A 939140 


85 2 


csr 
Cceanography Center 
946 — 


Chairman, Océanocraphy Department 
U.S. Naval Academ 


Anapolis, 


M 


D 2140 


Chief of Naval Research 


800 N. Qui 
Arlington; 


Chief of Naval 
Navai Ocean Research and Development 


ACTIVITY 


n 


VA 22217 


MEN EE Ste con 


Bay St. LOUIS, 


Commander 


Research 


et 
L 


40539522 


(Code 380) 


Oceanographic Systems Pacific 


Box 1390 


Pearl Harbor, 


Commanding Office 


Naval Eastern Deren aphy 


Naval Air 


37 


ation 


Norfolk,VA 23511 


Meee Larry 


Department of Oceanographv 


Naval Post 
Monterey, 


reaker, 


raduate 
A 93940 


n1 795560 


Cede 68 


School 


136 


Center 














Thesis 
N36135 Nelson 
21 Estimation and map- . 


ping of cloud and 
rainfall areas with 
an interactive 
computer. 





